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ABSTRACT 

The  cascade  algorithm  that  is  used  for  extended  surface  analysis  depends  on  a 
new  parameterization  called  the  thermal  transmission  matrix  to  represent  a  single  fin. 
This  thermal  transmission  matrix,  which  is  intended  to  replace  the  more  familiar  fin 
efficiency  as  a  design  and  analysis  parameterization,  is  a  linear  transformation  that 
maps  conditions  of  heat  flow  and  temperature  at  the  fin  tip  to  heat  flow  and 
temperature  conditions  at  the  fin  base.    The  cascade  algorithm  was  derived  by 
resorting  to  an  analogy  between  a  fin  and  the  electrical  transmission  line.   The 
cascade  algorithm  permits  a  fin  to  be  subdivided  into  many  subfins  each  having  a 
thermal  transmission  matrix  and  then  the  individual  transmission  matrices  for  each  of 
the  subfins  can  be  used,  via  a  simple  matrix  product  to  form  an  overall  equivalent 
thermal  transmission  matrix  for  the  entire  fin.    This  thesis  develops  a  thermal 
transmission  matrix  for  the  radiating  rectangular,  trapezoidal,  and  triangular  fins  both 
for  the  free  space  and  non-free  space  environments.   Test  cases  have  been  run  and 
their  solutions  exactly  match  those  contained  in  the  literature.   The  thesis  concludes 
with  optimization  studies  for  each  profile  considered  where  it  is  observed  that  simple 
algebraic  equations  can  be  employed  to  describe  the  optimum  geometry. 
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I.    INTRODUCTION 

Since  the  advent  of  the  space  age,  extended  surfaces  have  been  extensively 
studied  as  a  means  for  heat  rejection  from  all  sorts  of  vehicles.    Because  weight 
optimization  is  of  the  utmost  importance  in  spacecraft  design,  it  is  desirable  to  have 
an  accurate  technique  that  will  rapidly  optimize  the  geometry  of  the  extended  surface 
used  for  rejecting  the  heat  absorbed  on  or  generated  within  the  vehicle.    To  this  point 
in  time,  the  optimization  problem  involved  the  time  consuming  solutions  of  second 
order  non-linear  differential  equations,  by  numerical  methods  or  by  cumbersome  and 
possibly  inaccurate  graphical  methods.    A  more  advantageous  method  for  determining 
the  minimum  mass,  maximum  radiant  heat  rejection  extended  surface  geometry  has 
been  needed. 

The  purpose  of  this  study  is  to  employ  the  new  parameterization  called  the 
thermal  transmission  matrix  to  represent  a  single  fin  to  evaluate  the  performance  of 
several  longitudinal  fins.   Each  fin  studied  is  subdivided  into  100  subfins  and  the 
thermal  transmission  matrix  is  calculated  for  each  subfin.   Then,  by  a  simple  matrix 
multiplication  of  the  individual  subfin  thermal  transmission  matrices,  the  overall 
transmission  matrix  is  computed.    Not  only  will  this  method  produce  the  base  and  tip 
conditions  of  the  overall  fin  entity,  but  the  temperature  and  heat  flow  at  any  point 
between  the  fin  base  and  fin  tip  can  be  determined.    This  iterative  method  considers 

13 


the  spacecraft  radiant  heat  rejection  problem  in  both  the  free  and  non-free  space 
environments.    Optimization  must  be  considered  and  simple  algebraic  equations  that 
describe  the  fin  geometry  are  formulated  for  the  longitudinal  fins  of  rectangular, 
trapezoidal,  and  triangular  profiles. 

This  new  approach  to  fin  synthesis  and  analysis  permits  an  exact  solution  to  any 
degree  of  accuracy  that  the  user  wishes  to  specify  and  produces  the  result  in  a 
computational  efficient  manner. 
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H.    THEORETICAL  PRESENTATION 

A.      GENERALIZED  DIFFERENTIAL  EQUATION 

There  are  three  distinct  modes  by  which  heat  flows  from  one  point  to  another. 
A  rather  simplistic  definition  of  each  can  be  stated.    Conduction  is  the  mode  by  which 
heat  flows  through  a  material  by  molecular  motion.   If  the  heat  flows  from  a 
confining  surface  to  a  moving  fluid,  the  mode  is  known  as  convection.    If  no 
intervening  medium  is  present  and  the  thermal  energy  is  transported  from  one  point  to 
another  in  the  form  of  electromagnetic  waves,  the  mode  is  known  as  radiation  or 
radiant  heat  transfer. 

The  flow  of  heat  by  conduction  is  proportional  to  the  temperature  gradient  and 
cross  sectional  area  perpendicular  to  the  direction  of  heat  flow. 

At  ,  dT\x)  (2.1) 

Insertion  of  a  proportionality  constant  yields 

qcd--k  A(x)^&  (2.2) 

ax 

where  the  minus  sign  assures  a  positive  heat  flow  in  the  presence  of  a  temperature 
gradient  which  must  be  negative.   This  serves  to  define  the  thermal  conductivity  of 
the  material  written  as  k(T)  because  it  is  a  material  property  almost  always  a  function 
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-  Aix)  =^9- 
dx 

of  temperature.    Thus,  the  heat  flow  at  x  normal  to  a  surface  of  area  A(x)  at  a 
temperature  T(x)  is  written  as 

qcd  =  ~  k(T)  A(x)  —f-£ 
ax 

In  convection  the  heat  flow  from  a  confining  surface  at  temperature  T(x)  to  a 
fluid  whose  bulk  temperature  is  Ts  will  be  proportional  to  the  temperature  difference 
and  the  amount  of  surface  area. 

qcv  «  S(x)  [  nx)  -  Ts  ]  (2.5) 

Insertion  of  a  proportionality  constant  permits  the  representation 

qcv  =  h^  S(x)  [  T(x)  -  Ts  }  (2.6) 

and  this  serves  to  define  the  convective  heat  transfer  coefficient  hcv  as 

h      =  SSL (2.7) 

CT       S(x)  [  1\x)  -  Ts  ] 
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and  thus  Newton's  Law  of  Cooling  is  written  as 

<?cv  =  K^  5W  [  nx)  ~  Ts]  (2.8) 

where  the  convective  heat  transfer  coefficient  hcv  is  written  as  a  function  of  position 
and  temperature. 

Radiation  requires  no  intervening  medium  and,  for  a  perfect  emitter  and 
absorber  in  full  view  of  one  another,  the  familiar  Stefan-Boltzmann  law  shows  that  the 
heat  flow  is  proportional  to  the  surface  area  and  the  difference  of  the  fourth  powers  of 
the  absolute  temperature. 

qr  =  a  S(x)  [  T*  -  T24  ]  (2.9) 

The  proportionality  constant  used  here  is  the  Stefan-Boltzmann  constant, 
a  =  5.66961  x  10"8  watt  /  m2  K4. 

If  the  surfaces  are  not  perfect  emitters  and  not  perfect  absorbers  and  they  are 
not  in  full  view  of  one  another,  the  departure  from  perfection  is  via  the  emissivity 
factor  F£  and  the  shape  or  arrangement  factor  FA. 

qr  =  a  FAF9S(x)  [  7\4  -  T2*  ]  (2-10) 

The  problem  of  a  finned  surface  radiating  to  outer  space  simplifies  the  foregoing 
equation  (2.10).   The  temperature  of  the  fin  becomes  T(x)  and,  because  outer  space  is 
at  about  2  K,  T2  can  be  assumed  to  be  equal  to  0.    [Compare  T(x)4  =  2004  to  Ts4  = 
24].    Moreover,  FA  is  equal  to  1,  the  value  of  a  small  body  in  a  large  enclosure  (small 
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Figure  2.1   Energy  Balance 

fin  in  a  large  enclosure,  outer  space).    Finally  the  emissivity  factor  for  this  case 
(small  body  in  a  large  enclosure)  is  F,  =  e(x,T).   The  emissivity  e  is  a  property  of  the 
radiating  fin  surface  and  temperature.    Thus  equation  (2.10)  is  rewritten  as 


qr  =  a  e(x,T)  S(x)  T(xY 


(2.11) 


The  difference  between  the  rates  of  heat  entering  and  leaving  the  element  in 
Figure  2.1  by  conduction  (dqcdb  -  dqcda)  must,  in  the  steady  state,  be  equal  to  the 
rates  of  heat  dissipated  from  the  face,  L  dx  by  convection  (dqcvl  and  dqcv2)  and 
radiation  (dqrl  and  dqr2)  plus  any  rates  of  heat  that  arrive  from  external  sources  (dq^ 


18 


and  dq^). 

The  presence  of  two  terms  for  dqcv,  dqr,  and  dq,  permit  unequal  values  on  the 
two  faces.    Thus  with 


d 

dx 


d^cdA  ~  ^cd,b  z   — 


the  two  convective  dissipation  terms  are 


k(x)  b(x)  L 


dT 
dx 


dx 


^cvl    +   *la,2    =   KMV   L  &  [   T  ~   TSI   } 


+  hjtf)  Ldx[T-TS2] 


and  the  two  radiation  dissipations  terms  with  FA  =  1  and  ¥(  =  e(x,T)  are 


(2.12) 


(2.13) 


^rl    +    d(lr2    =    °    e(X'7)    L   d* 


4       t4 


r4  -z 


SI 


+  a  e(x,T)  L  dx[T4  -Ts\ 
The  energy  balance  with  TS1  =  Ts2  =  0  is 


(2.14) 


d_ 
dx 


k(x)  6(x)  L  — 
dx 


dx  =  [  h^ixj)  +  hjixj)  ]LdxT 


+  o  e(x,T)  Ldx  T4  +  qa  +  $r  (2.15) 


'.2 


In  space,  if  either  face  of  the  fin  is  subjected  to  an  external  heat  input  of  the 


form 


qi  =  a  L  E  dx 


(2.16) 
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where  E  consists  of  external  heat  inputs  from  solar,  Es,  and/or  terrestrial,  Er,  sources 
then  equation  (2.15),  the  steady  state  thermal  energy  balance  for  a  differential  element 
dx,  is  rewritten  as 


d_ 
dx 


k(x)  hix)  L  — 
dx 


d*  =  [  KM*1)  +  *«*W)  }LTdx 
+  o  e(x,T)  LT*dx+2aLEdx         (217) 


B.      LIMITING  ASSUMPTIONS 

Further  assumptions  must  be  applied  to  the  generalized  differential  equation 
(2.17)  before  a  mathematical  solution  can  be  attempted.   The  assumptions  employed 
by  Murray1  and  Gardner  are  used  for  a  starting  point  for  further  analysis.    The 
Murray-Gardner  Assumptions  listed  in  Table  1  are  essential  in  simplifying  the 
formulation  of  the  generalized  differential  equation  (2.17)  so  that  mathematical 
analysis  is  possible.    Consider  an  arbitrary  symmetrical  longitudinal  fin  shown  in 
Figure  2.2  of  length  (b-a),  depth  L,  and  width  5(x)  dissipating  radiant  heat  to  a 
constant  temperature  Ts  surrounding  environment. 

The  profile  area  Ap  is  defined  as 

x=b 

(2.18) 


Ap  =    /  6(x)  dx 


x=b-a 
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The  initial  conditions  at  the  fin  tip  are 

Ta  =  T(x=a)  (2.19) 

and 

qa  =  q(  x=a  )  (2.20) 

Applying  the  tenth  Murray-Gardner  assumption  (Table  1),  that  no  heat  leaves  from 
the  tip  of  the  fin,  equation  (2.20)  leads  to 


dn  x=a  >   =  0  (2.21) 


dx 

Two  more  boundary  conditions  are  taken  at  the  fin  base 


Tl  x=b  )  =  T.  (2.22) 


and 


q{  x=b  )  =  qb  (2.23) 

Applying  the  fifth  and  ninth  Murray-Gardner  assumptions  (Table  1),  allows  one  to 
ignore  the  radiant  heat  flow  from  the  edges.   The  surface  area  S(x)  perpendicular  to 
the  radiant  heat  flow  is  given  by 

S(x)  =  2  L  y  (2-24) 

The  face  area  S(x)  through  which  the  heat  is  dissipated  depends  on  the  arc  length  7,  a 
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function  of  the  fin  length  b. 

The  cross  sectional  area  A(x)  perpendicular  to  the  heat  flow  by  conduction  is 


A(x)  =  L  5(jc) 


(2.25) 


x-b 


Tb'6b  "  "  a 

.T(x),6(x)fq(x) 


A(x 


base,Tb,6  b,q 


Figure  2.2   Arbitrary  Symmetrical  Longitudinal  Fin3 

Because  the  second  Murray-Gardner  assumption  (Table  1)  states  that  the  material  is 
homogeneous,  both  the  emissivity  and  thermal  conductivity  are  assumed  to  be 
constant  and  uniform 
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TABLE  1 
THE  MURRAY-GARDNER  ASSUMPTIONS4 


1.  The  heat  flow  into  the  fin  and  the  temperature  at  any  point  on  the  fin 
remains  constant  with  time. 

2.  The  fin  material  is  homogeneous,  its  thermal  conductivity  is  the  same  in 
all  directions  and  remains  constant. 

3.  The  heat  transfer  coefficient  between  the  fin  and  the  surrounding  medium 
is  uniform  and  constant  over  the  entire  surface  of  the  fin. 

4.  The  temperature  of  the  medium  surrounding  the  fin  is  uniform. 

5.  The  fin  width  is  so  small  compared  with  its  height  that  temperature 
gradients  across  the  fin  width  may  be  neglected. 

6.  The  temperature  at  the  base  of  the  fin  is  uniform. 

7.  There  are  no  heat  sources  within  the  fin  itself. 

8.  Heat  transfer  to  or  from  the  fin  is  proportional  to  the  temperature  excess 
between  the  fin  and  the  surrounding  medium. 

9.  There  is  no  contact  resistance  between  fins  in  the  configuration  or  between 
the  fin  at  the  base  of  the  configuration  and  the  prime  surface. 

10.  The  heat  transferred  through  the  outermost  edge  of  the  fin  (the  fin  tip) 
is  negligible  compared  to  that  through  the  lateral  surfaces  (faces)  of  the  fin. 
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e(x,T)  =  e  (2-26) 

k(T)  =  k  (2-27) 

The  contribution  from  convective  dissipation  from  the  faces  of  the  fin  is  ignored 
because  the  arbitrary  symmetrical  longitudinal  fin  is  operating  in  a  space  environment. 

C.      GENERALIZED  DIFFERENTIAL  EQUATION  FOR  A  FIN 

Combining  the  generalized  differential  equation  (2.17)  and  equations  (2.24) 
through  (2.27),  and  assuming  that  there  are  no  external  heat  inputs  (^  =  qa  =  0), 
the  steady  state  thermal  energy  balance  for  the  differential  element  shown  in  Figure 
2.2  is 

A(k  I  *(A*M)  -  2  L  a  e  T(x)4  *L  =  0  (2.28) 

dx{  dx   )  dx 

The  temperature  excess  0(x)  is  defined  as  the  temperature  difference  between  the 
temperature  at  point  x  and  the  constant  temperature  of  the  surrounding 
environment  Ts. 

0(x)  =  T(x)  ~  Ts  (2.29) 

and  from  this,  it  is  seen  that 

dT{x)  -  dQ(x)  (2-30) 

Applying  the  temperature  excess  equations  (2.29)  and  (2.30),  the  equation  for  the 
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steady  state  thermal  energy  balance  for  the  differential  element  is  rewritten  as 


d_ 
dx 


&(*) 


dd(x) 
dx 


2  hr  dv 

-  6(x)  =2- 

k  dx 


(2.31) 


where  hr  is  defined  as  the  radiative  heat  transfer  coefficient.    Because  Ts  =  0, 


hr  =  o  e  Q(x)2 


(2.32) 


A  fundamental  concept  of  calculus  is  that  the  length  of  arc  d7  in  a  cartesian 
coordinate  framework  is 


dy  -  [  dx1  +  dy2}  m 


(2.33) 


The  derivative  of  the  arc  length  is 


dy 

dx       Ax-o  Ax 


im*X  -± 


N 


'  •  '21 


(2.34) 


The  +  or  -  sign  is  to  be  taken  according  as  the  arc  length  7  increases  or  decreases  as 
x  increases. 

Combining  the  equation  for  the  steady  state  thermal  energy  balance  differential 
element  (2.31)  and  the  arc  length  equation  (2.34),  the  generalized  differential  equation 
for  a  longitudinal  fin  of  arbitrary  symmetrical  profile  is 


d_ 
dx 


Hx) 


dQ(x) 
dx 


2  h 


T-  eoo  [1  +  t(X)2} m  =  0 


(2.35) 


Equation  (2.35)  is  a  second  order,  linear,  homogeneous  differential  equation  with  the 
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boundary  conditions 

dQ(  x=a  >  =  0  (2.36) 

dx 

and 

8(  x=b  )  =  db  (2.37) 

D.      LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE 

The  longitudinal  fin  of  rectangular  profile  is  illustrated  in  Figure  2.3.   It  should 
be  noted  that  the  origin  of  the  x  axis  is  located  at  the  tip  of  the  fin  and  not  at  the  base 
of  the  fin.   Because  the  width  is  constant  over  the  entire  profile,  the  width  function 
<5(x)  is  simply 

fi(x)  =  6,  =  8fl  (2.38) 

The  profile  area  Ap  is  written  as 

x=b 


Ap  =   fb(x)  dx  =  bbb  (2-39) 


*=o 


Inserting  the  width  function  (2.38)  into  the  generalized  differential  equation  for  a  fin 
equation  (2.35)  yields 


s/eoo  -  i*t  eoc)  =  o  (2.40) 

b  dx2  k 
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or 


d2Q(x) 


where 


-  m2  0(x)  =  0  (2.41) 


m  = 


2  K  (2.42) 


\  k  b 


b 


The  radiative  heat  transfer  coefficient  hr  is  assumed  to  be  a  constant  for  a  very  small 
subfin  where  Tb  «  Ta. 

The  general  solution  of  the  generalized  differential  equation  (2.41)  is 

0(jO  =  C,  emx  +  C2  e-'"x  (2-43) 

where  Q  and  C2  are  arbitrary  integration  constants  that  are  evaluated  by  applying  the 
initial  condition  of  equation  (2.22). 

8(6)  -  %  -  C1  em  b  +  C2  e-m  b  (2-44) 

The  heat  flow,  qb,  into  the  base  of  the  fin  is  found  by  inserting  the  general 
solution  of  the  generalized  differential  equation  (2.43)  and  the  cross  sectional  area 
perpendicular  to  the  heat  flow  by  conduction  (2.25)  into  Fourier's  Law  of 
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Figure  2.3   Longitudinal  Fin  of  Rectangular  Profile 

Conduction,  equation  (2.4). 


q(x)  =  Lbbk(  j-[Cl  emx  +  C,«""x]' 


(2.45) 


Simplifying  and  evaluating  at  x  =  b  with  the  initial  condition  of  equation  (2.23), 
equation  (2.45)  becomes 


qb  =  L  hbkm  [Cxemb  -  C2  e~m  b  ) 


(2.46) 
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Solving  for  C,  and  C2  yields 


„-mb                  0-mb 
Cl    =    " Qb    +    ~ 

1  2        b       2  YL 


% 


(2.47) 


and 


,m  b 


,m  b 


2  yt 


<ib 


(2.48) 


where 


yA  =  I  (  2  ftf  it  6,)  ^  (2.49) 

is  defined  as  the  characteristic  thermal  admittance  Yb  of  the  fin.    Likewise  the 
characteristic  thermal  impedance,  Zb,  of  the  fin  is  the  reciprocal  of  Yb 


z*  ■  j  -  { ( 2 »,  * », ) 


-1/2 


(2.50) 


An  examination  of  equations  (2.43)  through  (2.48)  reveals  that  only  two  of  the 
four  variables  are  independent  and  that  the  specification  of  any  two  variables  allows 
for  the  determination  of  the  remaining  two.    Combining  equations  (2.43)  through 
(2.48)  yields  the  two  port  inverse  thermal  transmission  matrix  for  a  longitudinal  fin  of 
rectangular  profile 


6 


cosh  nib      -Zb  sinh  mb 


-Yb  sinh  mb 


cosh  mb 


(2.51) 
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where  the  inverse  thermal  transmission  matrix  is  defined  as 


cosh  mb      -Zb  sinh  mb 


-Yb  sinh  mb 


(2.52) 


cosh  mb 

The  thermal  transmission  matrix  is  defined  as  the  inverse  of  the  inverse  thermal 
transmission  matrix. 


T  =  T1  = 


A      B 
C     D 


(2.53) 


where  A,  B,  C,  and  D  are  the  elements  of  the  thermal  transmission  matrix.    The 
thermal  transmission  matrix  representation  is 


&; 

cosh  mb 

Zb  sinh  mb 

X 

qt>. 

Yb  sinh  mb 

cosh  mb 

s. 

(2.54) 


E.      LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE 

The  longitudinal  fin  of  trapezoidal  profile  is  shown  in  Figure  2.4.    Note  that  the 
origin  of  the  axis  is  not  at  the  tip  of  fin  but  located  to  the  right  of  the  tip 
where  x  =  0. 


The  width  function  <5(x)  for  this  profile  is 


b(x)  =  2  x  tan  4) 


(2.55) 
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where  <f>  is  defined  as 


4>  =  arctan 


2  (b  -  a) 


(2.56) 


The  profile  area  A   is 


x=b 

Ap=  fmdx=V^(6b  +  6a) 


(2.57) 


Inserting  the  width  function  (2.55)  into  the  generalized  differential  equation  for  a  fin 
[Equation  (2.35)]  yields 


2  x  tan«|.)  *§&  ♦  2  tanOM  Q®.  -  -1^-  6(x)  =  0        (2.58) 
dx1  dx  k  cos((J>) 


or 


d2Q(x)       dQ(x)  , 


dx1  dx 


+  -^^  -  mL  Q(x)  =  0 


(2.59) 


where 


m  = 


k  sin(4>) 


1/2 


(2.60) 


Equation  (2.59)  represents  a  second  order,  linear,  homogenous  differential  equation. 
The  general  solution  is  expressed  in  terms  of  modified  Bessel  functions.   The  general 
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Figure  2.4  Longitudinal  Fin  of  Trapezoidal  Profile5 

solution  to  equation  (2.59)  is 

Q(x)  =  Cj  Io{  2  m  x  1/2 )  +  C2Ko(2  m  x112)  (2.61) 

The  following  transformation  is  made 


u  =  2  m  x  1/2 


(2.62) 


such  that  at  x  =  b 


ub  =  2  m  b  1/2 


(2.63) 
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and  at  x  =  a 


u.  =  2  m  a  m 


(2.64) 


After  applying  the  boundary  conditions  of  equations  (2.36)  and  (2.37)  and 
making  the  transformations  of  equations  (2.63)  and  (2.64),  the  general  solutions 
become 


e*  -  c,  W  +  c2  K0(Ub) 


(2.65) 


and 


<lb 


2  k  bbL  m2 


[C,  W  -  C2  *,(«,)] 


(2.66) 


Simultaneous  solution  of  the  general  solution  equations  (2.65)  and  (2.66)  yield 


C,  =«* 


*W  6, 


2  k  bbL  m2 


*W  <l» 


(2.67) 


and 


C2  =«* 


W  ©I 


2khbLm2 


W  <?> 


(2.68) 


Substituting  and  rearranging  the  general  solution  equations  (2.65)  and  (2.66)  and  the 
equations  for  Q  and  C2  (2.67)  and  (2.68),  the  elements  of  the  thermal  transmission 
matrix  for  the  longitudinal  fin  of  trapezoidal  profile  are  written  as: 
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A  =  ua  [  J^  Ke(uJ  *  I0(Ub)  Kfiij] 


(2.69) 


B  = 


2  k  8, 


— -    [WW  -  /0(ua)  K0(Ub)] 
L  m 


(2.70) 


C  =  2  k  bLm2 


fu\ 


\U*J 


[  Ifr}  Kx{ua)  -  I^UJ  K^ujl  (2.71) 


0  =  H  i  W  *W  +  AW  W] 


(2.72) 


F.      LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE 

The  longitudinal  fin  of  triangular  profile  shown  in  Figure  2.5  is  a  special  case  of 
the  longitudinal  fin  of  trapezoidal  profile.   There  are  many  similarities  between  the 
longitudinal  fin  of  triangular  and  trapezoidal  profiles.    The  differential 
equation  for  the  longitudinal  fin  of  triangular  profile  is  the  same  as  the  differential 
equation  for  the  longitudinal  fin  of  trapezoidal  profile  equation  (2.59). 


d2Q(x)       dQ(x) 


+  -  m2  Q(x)  =  0 

dx2  dx 


(2.73) 
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Figure  2.5   Longitudinal  Fin  of  Triangular  Profile6 


where 


and 


1/2 


m  = 


k  sin(<J>) 


(2.74) 


(J)  =  arctan 


2  b 


(2.75) 
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The  profile  area  A   is 


x=b 

b 


Ap  =  fb(x)dx  =  |  6,  (2-76) 


x=0 

However  the  evaluation  of  the  arbitrary  constants  from  the  boundary  conditions 
differ  between  the  triangular  and  the  trapezoidal  profiles.    As  the  value  of  x 
approaches  zero  the  temperature  excess  at  the  tip  6&  must  be  equal  to  some  finite 
value,  therefore 

lo(  u=0  )  =  1  (2.77) 

In  order  to  physically  obtain  a  finite  temperature  at  xa  =  ua  =  0,  C2  in  equation 
(2.61)  must  be  equal  to  zero  because  K2(0)  is  unbounded.   The  unbounded  condition 
allows  cancellation  of  the  K0(u)  terms  from  the  general  solutions  for  temperature 
excess  and  heat  flow. 

After  cancellation  and  rearranging  the  general  solutions  for  the  longitudinal  fin 
of  trapezoidal  profile  equations  (2.65)  and  (2.66),  the  general  solution  for  the 
longitudinal  fin  of  triangular  profile  becomes 


=  2  k  bb  L  m2   W  e  (2.78) 


As  a  result  of  the  zero  tip  width,  the  edge  parameters  do  not  depend  on  the  base 
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parameters.    The  thermal  input  admittance  element  for  the  longitudinal  fin  of 
triangular  profile  is 


y   =    q>   =   2  k  bbLm2    Ix{uJ)  (2J9) 


Qb  ub  'oW 


m.    CASCADE  ANALYSIS  ALGORITHM 

A.      GENERALIZED  TWO  PORT  TRANSMISSION  LINE  MODEL 

The  radiating  fin  is  a  distributed  heat  flow  configuration  which  is  analogous  to 
the  transmission  line  which  is  a  distributed  electrical  network.    The  general  two  port 
transmission  line  model  is  shown  in  Figure  3.1.   The  model  is  described  by  two  input 
terminals.    At  point  "a"  there  is  an  applied  voltage  Va  and  a  current  Ia  flows  into  the 


Transmission 
Line 

'■     v. 

vb       '> 

Figure  3.1   Two  Port  Transmission  Line  Model 

model.   There  is  also  a  pair  of  terminals  at  point  "b"  where  a  voltage  Vb  occurs  and  a 
current  Ib  flows  out  of  the  terminals.    The  two  port  model  of  the  transmission  line  can 
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consist  of  any  arrangement  of  passive  elements  within  the  two  port  but  with  no  active 
current  or  voltage  sources  within  the  two  port. 

It  can  be  shown  from  network  theory  that  the  two  port  transmission  line  model 
in  Figure  3.1  can  be  represented  by  the  following  equations: 


Vm-AVt*BIb 


(3.1) 


Ja  =  C  Vb  +  D  Ib 


where  the  determinant  of  the  ABCD  matrix  defined  by 


(3.2) 


= 

A  B 
C  D 

is  equal  to  unity. 


A  D  -  B  C  =  1 


(3.3) 


(3.4) 


The  values  A,  B,  C,  and  D  are  constants  and  are  defined  as  the  elements  of  the 
transmission  matrix. 

The  electrical  transmission  line  analogy  can  be  directly  applied  to  an  equivalent 
thermal  transmission  line  model  where  the  heat  flow  is  analogous  to  current  and  the 
temperature  is  analogous  to  voltage.   Equations  (3.1)  and  (3.2)  are  rewritten  in  the 
equivalent  thermal  transmission  line  model  as: 


Ta  =  A  Tb  *  B  qb 


(3.5) 
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ga  =  C  Tb  +  D  qb 


(3.6) 


B.      CASCADE  ALGORITHM 

The  thermal  transmission  parameters  are  very  useful  in  describing  two  port 
networks  connected  in  a  cascaded  arrangement.  It  is  an  exercise  in  elementary  two 
port  theory.    Examine  the  two  equations  for  the  two  networks  shown  in  Figure  3.2. 


= 

A,       B, 

V 

(3.7) 


and 


V 

.q\ 

= 

\          B2 

C2      D2 

V 

(3.8) 


Continuity  for  the  configuration  in  Figure  3.2  dictates  that 


(3.9) 


and 


qb2    =    ^ 


(3.10) 


which  can  be  represented  as 


(3.11) 


The  overall  thermal  transmission  parameter  matrix  for  the  two  networks 
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cascaded  in  Figure  3.2  is  by  definition 


A 

B 

\ 

Br" 

\ 

B2 

C 

D 

Px 

DL 

C2 

D2. 

(3.12) 
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Figure  3.2   Cascaded  Network 

and  the  thermal  transmission  matrix  for  a  cascaded  network  of  n  networks  is 


A 

B 

\ 

K 

A2 

*2 

'*n 

Bn 

C 

D 

<k 

Di. 

C2 

D2. 

cn 

Dn. 

(3.13) 


Equation  (3.13)  is  an  expression  of  what  is  referred  to  as  the  cascade  algorithm. 


C.      COMPUTER  ANALYSIS 

The  cascade  algorithm  expression  (3.13)  can  be  applied  to  a  fin  divided  into  n 
sub  fins  each  having  the  same  length.    Converting  the  temperature  excess  notation 


back  to  the  temperature  notation  (Ts  *  0),  the  cascade  algorithm  expression  (3.13)  is 


written  as 


= 

\      Bi 

Cx      D± 

A2       B2 

C2         D2. 

\        Bn 

q*. 

(3.14) 

The  base  temperature  Tb  and  the  base  heat  flow  qb  are  the  conditions  at  the  base  and 
the  tip  temperature  Ta  and  tip  heat  flow  qa  are  the  conditions  at  the  edge. 

A  problem  arises  because  the  tip  temperature  Ta  and  the  base  heat  flow  qb  are 
not  known  or  not  specified.    The  cascade  algorithm  expression  (3.14)  can  not  be 
solved  directly  because  it  is  essentially  a  linear  transformation  that  maps  tip  conditions 
to  base  conditions.    However,  a  solution  can  be  calculated  indirectly.   If  one  makes 
an  initial  estimation  for  the  tip  temperature  Ta,  it  can  be  used  to  calculate  the  radiant 
heat  flow  qa  out  of  the  tip  of  the  fin. 


qa  =  o  e  5a  L  Ta 


(3.15) 


Thus  Ta  and  qa  are  the  starting  values  for  the  computations  using  the  cascade 
algorithm  expression  (3.14).   At  the  conclusion  of  the  process,  a  comparison 
involving  all  n  subfins  is  then  made  between  the  known  or  specified  base  temperature 
Tb  and  the  calculated  base  temperature.   If  the  comparison  indicates  an  error,  the  tip 
temperature  Ta  is  changed  and  the  iterative  process  is  repeated  and  continues  until  the 
error  between  the  presumed  known  base  temperature  Tb  and  the  calculated  base 
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temperature  are  within  an  acceptable  range.    Figure  3.4  provides  a  complete  flow 
chart  of  the  cascade  algorithm. 

Further  amplification  concerning  the  radiative  heat  transfer  coefficient  needs  to 
be  made.    Recalling  equation  (2.32)  where 


hr  =  a  e  T{x) 


(3.16) 


it  is  not  clear  as  to  what  temperature,  T(x),  should  be  applied  to  the  computation  of  h, 
and  hence  to  the  cascade  algorithm  expression  (3.14).   The  left  and  right  faces  of  the 
differential  element  in  Figure  3.3  have  an  average  temperature. 


Ta    +    Tb 


(3.17) 
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Figure  3.3   Differential  Element 
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and  as  long  as  Ta  and  Tb  are  not  vastly  different,  equation  (3.16)  can  be  written  using 


hT  =  o  e  Tj  (3.18) 

The  environmental  parameter  V  is  introduced  for  the  case  in  which  there  are 
portions  of  the  surrounding  environment  at  some  temperature  other  than  Ts  «  0.    An 
example  of  this  is  in  the  interception  of  thermal  radiation  from  a  nearby  solar  panel  at 
some  temperature  not  near  absolute  zero.    The  environmental  parameter  is  defined  as 

V  =   1 (3.19) 

2   o  €  Tb* 

where  the  parameter  K2  allows  for  the  application  of  external  heat  inputs  from  solar, 
Es,  or  terrestrial,  ET,  sources  or  both. 

K2   =  E  a  (3.20) 

The  environmental  parameter  is  equal  to  zero  for  the  free  space  case  and  to  some 
number  greater  than  zero  but  less  than  one  for  the  non  free  space  case.    Substituting 
equation  (3.19)  into  (3.18),  the  radiant  heat  transfer  coefficient  becomes 


hz  =  o  e 


T    3  -       Kl 

3.V  *y       rp 


(3.21) 


External  heat  inputs  (q^  and  q^)  from  solar  input  Es  and  terrestrial  input  Ef  are  now 
applied  to  the  cascade  algorithm  expression  (3.14)  through  the  application  of  the 
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environmental  parameter  V. 

D.      PROGRAM  DESCRIPTION 

The  computer  code  consists  of  three  analysis  programs  for  each  of  the  profiles 
considered,  and  a  MAIN  Program  that  is  used  for  general  housekeeping.    The 
PERFORMANCE  Analysis  Program  evaluates  an  entered  set  of  parameters  returning 
the  characteristics  (volume,  profile  area  Ap,  profile  number  PN,  base  heat  flow  qb, 
edge  temperature  Ta,  and  the  fin  efficiency  rj)  of  the  fin  along  with  the  thermal 
transmission  parameter  matrix.    The  PROFILE  AREA  Analysis  Program  also 
evaluates  an  entered  set  of  parameters  but  it  returns  a  data  set  of  100  fins  having  a 
constant  profile  area  Ap  with  the  base  fin  width  5  and  length  b  which  are  used  as  the 
input  variables.    The  OPTIMIZE  Analysis  Program  is  similar  to  the 
PERFORMANCE  Analysis  Program  with  the  exception  that  it  returns  the  optimum 
fin  geometry  which  results  in  the  minimum  mass  and  maximum  base  heat  flow  fin. 

The  MAIN  Program  provides  menu  driven  access  to  the  various  analysis 
programs.    The  user  has  the  ability  to  set  the  computational  main  program  error 
coefficient.   The  user  also  has  the  ability  to  set  the  save  data  switch  which  controls 
data  flow  into  an  ASCII  formatted  disk  file.   Commas  are  used  for  the  data  delimiter 
which  allows  file  importation  into  MATLAB,  Word  Perfect,  Draw  Perfect,  and  Math 
CAD  for  additional  analysis  or  for  plotting.    An  additional  feature  provided  is  the  unit 
conversion  module.   Table  2  list  the  available  units  that  are  provided  to  the  user  by 
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the  MAIN  Program.    Figure  3.4  provides  a  flow  chart  of  the  cascade  algorithm 
expression  (3.14)  which  is  common  to  all  of  the  analysis  programs. 

The  PERFORMANCE  Analysis  Program  requires  the  following  input  conditions 
to  be  entered: 

•  base  temperature 

•  density 

•  emissivity 

•  environmental  parameter 

•  fin  base  width 

•  fin  edge  width1 

•  fin  depth 

•  fin  length 

•  thermal  conductivity 

The  program  calculates  and  returns  the  following  values: 

•  base  heat  flow 

•  edge  temperature 

•  fin  efficiency 

•  profile  area 

•  profile  number 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  fin  efficiency  rj  and  the  profile  number  PN  are  extensively  referred  to  in 
literature  and  are  provided  for  use  in  making  comparisons.    The  fin  efficiency  and  the 
profile  number  are  defined  as 


Tl    =  2S _  (3.22) 

2  o  e  b  L  Tb4 


'Applies  for  only  the  longitudinal  fin  trapezoidal  profile 
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TABLE  2 


UNITS  AVAILABLE  IN  THE  MAIN  PROGRAM 


Admittance 

•  watts  /  °C,   BTU  /  hr  °F 

Area 

•  m2,   cm2,   ft2,   in2,   mm2 

Heat 

•  watts,  kw,   BTU  /  sec,   BTU  /  hr 

Impedance 

•  °C  /  watts,    °F  hr  /  BTU 

Length  and  Depth 

•  m,   cm,   ft,   in,    mm 

Mass 

•  kg,   grams,   lbm,   oz 

Temperature 

•  K,    °C,    °F,    °R 

Thermal  Conductivity 

•  watts  /  m  K,   BTU  /  ft  hr  °F 

Width 

•  m,   cm,    ft,   in,   mm 

PN  =    2  q  e  _     T*      (trapezoidal)    (3.23) 
k  5. 


PN  = 


q  e  b2  Tb3    (rectangular)     (3.24) 


If  the  save  switch  is  switched  on,  the  following  differential  element  conditions  are 
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written  to  a  disk  file: 

•  base  heat  flow 

•  base  temperature 

•  distance  from  the  tip  of  the  fin 

The  OPTIMIZE  Analysis  Program  requires  the  following  input  conditions  to  be 
entered: 

•  base  temperature 

•  density 

•  emissivity 

•  environmental  parameter 

•  fin  edge  width2 

•  fin  width 

•  thermal  conductivity 

The  program  calculates  and  returns  the  following  optimum  values: 

•  base  heat  flow 

•  edge  temperature 

•  fin  base  width 

•  fin  efficiency 

•  fin  length 

•  profile  area 

•  profile  number 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  program  saves  the  same  quantities  that  the  PERFORMANCE  Analysis  Program 
saves  depending  on  the  condition  of  the  save  switch. 

The  PROFILE  AREA  Analysis  Program  requires  the  following  input  conditions 


2Applies  for  only  the  longitudinal  fin  trapezoidal  profile. 
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to  be  entered: 


•  base  temperature 

•  density 

•  emissivity 

•  environmental  parameter 

•  profile  area 

•  starting  fin  base  width  value 

•  ending  fin  base  width  value 


•  taper  ratio3 


•  fin  depth 

•  thermal  conductivity 

The  program  calculates  and  writes  to  a  disk  file  the  following  parameters  for  100  fins 
having  the  same  profile  area  A,,  within  the  specified  width  range: 

•  base  heat  flow 

•  base  width 

•  base  temperature 

•  edge  temperature 

•  fin  efficiency 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  analysis  programs  all  use  the  thermal  transmission  matrices  computed  in 
Chapter  2  in  the  cascade  algorithm  with  the  exception  of  the  longitudinal  fin  of 
triangular  profile.   The  triangular  and  trapezoidal  profile  analysis  are  very  similar. 
The  analysis  program  uses  the  thermal  admittance  equation  (2.79)  to  compute  the  base 
conditions  for  the  first  subfin  having  a  thickness  of  5a.   The  subfin  is  located  at  the  tip 
of  the  fin.   The  base  conditions  are  then  applied  and  used  as  the  tip  conditions  in  the 
thermal  transmission  matrix  for  the  trapezoidal  profile  (2.69)  through  (2.72).   The 


3Applies  for  only  longitudinal  fin  of  trapezoidal  profile. 
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cascade  algorithm  is  then  used  to  calculate  the  thermal  transmission  matrices  for  the 
remaining  99  subfins  in  order  to  finally  determine  the  overall  thermal  transmission 
matrix  and  the  base  conditions.   The  trapezoid  that  is  being  used  in  the  calculation  has 
an  almost  zero  tip  width  compared  to  the  base  width  such  that 


6b  =  100  6a 


(3.25) 


Borland  Turbo  C  +  +  was  used  to  write  the  computer  code.    The  computer  code 
listing  appears  in  Appendix  A.    Full  compiler  optimization  was  selected  in  order  to 
decrease  the  run  times.    Average  run  times  on  an  IBM  486/487  machine  appear  in 
Table  3. 

TABLE  3 
AVERAGE  PROGRAM  RUN  TIMES 


PERFORMANCE 

PROFILE  AREA 

OPTIMIZE 

(sec) 

(sec) 

(sec) 

Rectangular 

4 

97 

4 

Trapezoidal 

4 

600 

4 

Triangular 

5 

625 

5 

IV.    PROGRAM  ANALYSIS  RESULTS  AND  DISCUSSION 

A.      PROGRAM  ERROR  TOLERANCE 

In  order  to  allow  for  meaningful  interpretation,  an  understanding  of  the  accuracy 
of  the  results  produced  from  the  cascade  algorithm  is  vital.    A  test  to  show 
convergence  is  used  to  illustrate  the  approximate  relative  error. 

Table  4  provides  an  illustration  of  how  the  error  coefficient  in  the  MAIN 
Program  affects  the  accuracy  of  the  program  results.    Table  4  is  produced  by 

TABLE  4 
PROGRAM  ERROR  TOLERANCE 


Program 

Error 

Ta 

•  % 

qb 

«  % 

Coeff  % 

(K) 

error 

(watts) 

error 

1.000 

282.18 

0.0106 

926.28 

2.0918 

0.500 

282.24 

0.0106 

931.83 

1.5052 

0.100 

282.19 

0.0071 

946.36 

0.0307 

0.050 

282.22 

0.0035 

946.94 

0.0920 

0.010 

282.22 

0.0035 

946.10 

0.0032 

0.005 

282.21 

0.0000 

946.05 

0.0021 

0.001 

282.21 

0.0000 

946.07 

0.0000 

running  the  same  input  parameters  many  times  while  decreasing  the  error  coefficient 
in  the  MAIN  Program.    The  PERFORMANCE  Analysis  Program  for  the  longitudinal 
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fin  of  rectangular  profile  was  used.    The  Table  4  results  are  applicable  for  any  of  the 
analysis  programs  because  of  the  commonality  of  the  cascade  algorithm.    As  the  error 
coefficient  in  the  MAIN  Program  decreases,  the  output  value  converges.    The 
approximate  relative  error  et  is  defined  as 


zr  * 


True  Value  -  Approximate  Value 


100  (4.1) 


True  Value 

The  exact  true  value  is  unknown,  however  the  value  produced  when  the  error 
coefficient  in  the  MAIN  Program  is  equal  to  0.001  can  be  assumed  to  be 
approximately  equal  to  the  exact  true  value. 

Examining  the  results  illustrated  in  Table  4,  one  can  conclude  that  when  the 
error  coefficient  in  the  MAIN  Program  is  less  than  or  equal  to  0. 1  the  approximate 
relative  error  observed  in  the  output  is  generally  much  less  than  the  error  coefficient 
in  the  main  program. 

B.       GRAPHICAL  EXAMPLES 

The  purpose  of  this  section  is  to  describe  and  provide  examples  on  how  each  of 
the  three  analysis  programs  can  be  used.    The  computer  code  listings  are  provided  in 
Appendix  A. 

1.       Example  One:   Longitudinal  Fin  of  Rectangular  Profile  Analysis 

A  longitudinal  fin  of  rectangular  profile  is  0.5  meter  long,  1.0  centimeter 
thick,  and  1.0  meter  deep.    The  fin  is  made  of  aluminum  (k  =  209.4  watt  /  m  K, 
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p  =  2.6  kg  /  dm3,  and  e  =  0.85).   If  the  base  temperature  is  400  K,  what  are  the  tip, 
midpoint,  and  base  conditions  (error  <  0.1%)  both  in  the  free  space  and  non  free 
space  environments? 

Because  the  geometry  of  the  fin  is  already  known,  the  PERFORMANCE 
Analysis  Program  is  used  for  the  evaluation.   The  save  switch  in  the  MAIN  Program 
must  be  placed  in  the  save  position  in  order  to  later  retrieve  the  midpoint  conditions 
from  the  disk  file.   The  computer  output  display  for  the  free  space  environment 
(V=0)  appears  in  Table  5. 

Figures  4.1  and  4.2  provide  a  detailed  graphical  representation  of  the 
required  conditions  both  in  the  free  space  and  non  free  space  environments.   The  base 
conditions  are  read  from  where  the  length  is  equal  to  50  centimeters  and  the  tip 
conditions  are  read  from  where  the  length  is  equal  to  0  centimeter.    In  the  free  space 
environment  the  tip,  midpoint,  and  base  conditions  for  this  particular  fin  geometry  are 
listed  in  Table  6. 
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TABLE  5 
PERFORMANCE  ANALYSIS  OUTPUT  DISPLAY 


THE  LONGITUDINAL  FIN 

OF 

1  RECTANGULAR  PROFILE 

Density 

= 

2.60  kg  /  dnT3 

Fin  Depth 

= 

1.  00  m 

Total  Mass 

= 

13.00  kg 

Fin  Length 

= 

0.50  m 

Fin  Width 

= 

1.00  cm 

Emissivity 

= 

0.85 

Profile  Area 

= 

50.00  cm~2 

Profile  Number 

= 

0.607 

Base  Heat  Flow 

= 

7  3  2.85  watts 

Fin  Efficiency 

= 

0.594 

Base  Temperature 

= 

400.07  K 

Edge  Temperature 

= 

324.45  K 

Thermal  Conductivity 

= 

209.40  watts  /  m  K 

Environmental  Parameter 

= 

0.00 

Transmission  Parameter  Matrix 

A  =  1.2287 

B  =  0.2583  deg  C  /  watt 

B  =  2.2374  watt  /  deg  C 

D  =  1.2842 

TABLE  6 
FIN  CONDITIONS  IN  FREE  SPACE 


Location                                Temp  q 

(K)  (watts) 

Base                                     400.07  732.85 

Midpoint                                341.66  292.33 

Tip                                        324.45  10.68 
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Figure  4.1   Heat  Dissipation  From  A  Rectangular  Profile  Fin 
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Figure  4.2   Temperature  Distribution  For  A  Rectangular  Profile  Fin 
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Equation  (3.4)  is  used  below  to  verify  that  the  determinant  of  thermal 
transmission  matrix  listed  in  Table  5  is  indeed  equal  to  unity. 


det 


1.2287   0.2583 
2.2374   1.2842 


1.0000  <4-2) 


2.       Example  Two:   Longitudinal  Fin  of  Rectangular  Profile  Area  Analysis 

A  longitudinal  fin  of  rectangular  profile  1.0  meter  deep  is  made  of 
aluminum  (k  =  209.4  watt  /  m  K,  p  =  2.6  kg  /  dm3,  and  e  =  0.85).   It  is  required 
to  remove  a  minimum  of  500  watts  of  heat  from  the  base  of  a  fin  at  a  temperature  of 
400  K.    Determine  the  geometry  of  the  fin  using  several  values  for  the  profile  area 

V 

The  PROFILE  AREA  Analysis  Program  is  used  because  the  geometry  of 

the  fin  is  unknown.    Figure  4.3  provides  a  graphical  representation  which  shows  the 

constant  profile  area  Ap  relationship  as  a  function  of  width  versus  heat  flow.    Any  of 

the  dimensions  listed  in  Table  7  would  meet  the  minimum  requirement  to  remove  the 

500  watts  of  heat  from  the  base  of  the  fin. 

Examining  Figure  4.3,  one  can  see  that  the  maximum  heat  flow  occurs  at 

a  particular  width  value.    More  heat  is  dissipated  as  the  profile  area  Ap  is  increased. 

The  value  of  the  width  at  the  point  of  maximum  heat  flow  also  increases  as 

the  profile  area  Ap  increases.    The  point  of  maximum  heat  flow  is  considered  as  the 

optimum  geometry  (maximum  heat  flow,  minimum  mass)  for  the  particular  profile 

area  Ap.   A  summary  of  the  optimum  geometry  points  extrapolated  from  the  data  files 
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FIGURE  4.3  Profile  Area  Analysis 
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TABLE  7 


PROFILE  AREA  ANALYSIS  POINTS 


\ 

5b 

b 

qb 

(cm2) 

(cm) 

(m) 

(watts) 

25 

0.334 

0.748 

507.4 

100 

0.490 

2.041 

633.7 

200 

0.724 

2.762 

771.2 

300 

0.334 

8.982 

523.1 

produced  from  the  PROFILE  AREA  Analysis  Program  appear  in  Table  8. 

TABLE  8 
OPTIMUM  GEOMETRY  POINTS  FOR  THE  RECTANGULAR  PROFILE 


AP 

^bo 

b0 

Mass0 

qbo 

qbo/m 

(cm2) 

(cm) 

(m) 

(kg) 

(watts) 

(w/kg) 

25 

0.646 

0.387 

6.50 

579.76 

89.19 

100 

1.543 

0.648 

26.00 

924.53 

35.56 

200 

2.440 

0.819 

52.00 

1165.83 

22.42 

300 

3.103 

0.967 

80.00 

1334.40 

16.68 

3.       Example  Three:    Comparison  of  Various  Profiles 

A  longitudinal  fin  1.0  meter  deep  having  a  profile  area  Ap  equal  to  50  cm2 
is  made  of  aluminum  (k  =  209.4  watt  /  m  K,  p  =  2.6  kg  /  dm3,  and 
e  =  0.85).   The  base  temperature  is  held  at  400  K.    Compare  the  relative 
performance  between  the  longitudinal  fins  of  rectangular,  trapezoidal,  and  triangular  profiles. 
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The  same  approach  is  used  as  in  the  previous  example.    A  longitudinal  fin 
of  trapezoidal  profile  having  a  taper  ratio  TR  =  0.5  is  used.    The  taper  ratio  is 
defined  as 

8 

TR  =  -2  (4.3) 

h 

Figure  4.4  provides  a  graphical  representation  which  shows  the  relationship  between 

the  three  profiles  being  considered. 

The  curve  in  Figure  4.4  produced  by  the  longitudinal  fin  of  trapezoidal 
profile  is  midway  between  the  curves  produced  by  the  longitudinal  fins  of  rectangular 
and  triangular  profiles.    Of  more  importance  is  the  fact  that  at  the  location  of  the 
optimum  geometry  points,  the  longitudinal  fin  of  triangular  profile  produces  the 
largest  heat  flow  while  the  longitudinal  fin  of  rectangular  profile  produces  the  smallest 
heat  flow.   The  optimum  geometry  points  extrapolated  from  the  disk  file  generated 
from  the  PROFILE  AREA  Analysis  Program  appear  in  Table  9.    Not  only  is  there  an 
optimum  fin  geometry  for  a  particular  profile  area  Ap  but  there  is  also  an  optimum 
profile  shape  which  produces  maximum  heat  flow. 

4.       Example  Four:   Optimum  Geometry  Comparison 

A  longitudinal  fin  2.0  meters  deep  (e  =  0.95)  is  required  to  remove  at 
least  1000  watts  from  a  spacecraft  in  a  free  space  environment  (V=0).   The  base 
temperature  is  held  at  400  K.    Using  the  list  of  materials  appearing  in  Table  10 
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Figure  4.4   Comparisons  of  Various  Profiles 
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TABLE  9 


OPTIMUM  GEOMETRY  POINTS  FOR  THE  VARIOUS  PROFILES 


Taper 

5b 

b 

Mass 

qb 

qb/m 

Profile 

Ratio 

(cm) 

(m) 

(kg) 

(watts) 

(W/kg) 

RECT 

1.0 

0.958 

0.522 

13.00 

732.47 

56.343 

TRAP 

0.5 

1.270 

0.524 

13.00 

772.93 

59.456 

TRI 

0.0 

1.309 

0.763 

13.00 

826.78 

63.598 

conduct  a  heat  flow  per  unit  mass  comparison  using  the  longitudinal  fins  of 
rectangular  and  triangular  profiles. 

TABLE  10 
MATERIALS  USED  IN  THE  OPTIMUM  COMPARISON8 


Chemical 

P 

k 

Material 

Symbol 

kg/dm3 

w  /  m  K 

Aluminum 

Al 

2.6 

209.4 

Copper 

Cu 

8.8 

384.0 

Iron 

Fe 

7.2 

58.0 

Magnesium 

Mg 

1.7 

157.0 

Molybdenum 

Mo 

10.2 

145.0 

Silicon 

Si 

2.3 

83.0 

Silver 

Ag 

10.5 

407.0 

Zinc 

Zn 

6.8 

140.0 
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The  OPTIMIZE  Analysis  Program  is  used  for  the  analysis.    The 
parameters  were  entered  individually  into  the  OPTIMIZE  Analysis  Program  and  the 
results  are  listed  in  Table  1 1 . 

TABLE  11 
OPTIMUM  COMPARISON  RESULTS 


*bo 

Material   cm 

7 

CI 

Mass0 
kg 

Sbo 

watts 

qbo/m 
w  /  kg 

Ag  RECT 
Ag  TRI 

0.40 
0.54 

21.73 
15.09 

45.63 
31.69 

997.86 
995.45 

21.86 
31.41 

Al  RECT 
Al  TRI 

0.79 
1.05 

42.24 
29.33 

21.96 
15.25 

998.96 
995.47 

45.49 
65.27 

Cu  RECT 
Cu  TRI 

0.43 
0.57 

23.03 
15.99 

40.54 
28.15 

997.94 
995.45 

24.61 
35.36 

Fe  RECT 
Fe  TRI 

2.84 
3.79 

152.49 
105.89 

221.11 
153.54 

1003.99 
995.75 

4.54 
6.48 

Mg  RECT 
Mg  TRI 

1.05 
1.40 

56.33 
39.12 

19.60 
13.61 

999.74 
995.49 

51.01 
73.14 

Mo  RECT 
Mo  TRI 

1.14 
1.51 

61.00 
42.36 

124.43 
86.41 

1000.00 
995.54 

8.04 
11.52 

Si  RECT 
Si  TRI 

1.99 
2.65 

106.56 
73.99 

49.66 
34.48 

998.41 
995.59 

20.11 
28.87 

Zn  RECT 
Zn  TRI 

1.18 
1.57 

63.17 
43.87 

85.92 
59.66 

1000.12 
995.50 

11.64 
16.68 
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Figure  4.5  shows  a  bar  graph  representation  of  the  heat  flow  per  unit  mass 
relationship  using  the  various  materials  listed  in  Table  10.    It  is  evident  from  Figure 
4.5  that  the  optimized  longitudinal  fin  of  triangular  profile  made  of  magnesium  is  the 
most  efficient  material  considered  of  those  listed  in  Table  10,  with  aluminum  being  a 
close  second.    The  melting  temperature  of  both  materials  (=  900  K)  is  far  above  the 
base  temperature  (400  K)  stated  in  the  problem  statement.    Comparing  Tables  10  and 
11  materials  having  a  small  value  for  the  density  (p  <  3.0  kg  /  dm3)  and  a  large 
value  for  the  thermal  conductivity  (  k  >  150  watt  /  m  K)  produce  large  values  for 
heat  flow  per  unit  mass. 

C.      DERIVATION  OF  WORKING  OPTIMUM  GEOMETRY  EQUATIONS 

As  shown  in  Tables  8  and  9  optimum  geometry  points  can  be  determined  using 
the  PROFILE  AREA  Analysis  Program.    The  overall  success  of  finding  the 
dimensions  for  the  optimum  geometry  using  the  PROFILE  AREA  Analysis  Program 
has  several  limitations.   The  first  limitation  is  that  the  starting  and  ending  fin  base 
width  values  must  be  known.    And  it  can  be  argued  that  a  second  limitation  is  that 
excessive  analysis  program  run  times  which  might  exceed  ten  minutes. 

Relationships  for  the  optimum  geometry  for  the  longitudinal  fins  of  rectangular, 
trapezoidal,  and  triangular  profiles  would  remove  both  limitations.    Kern  and  Kraus9 
and  Chung  and  Nguyen10  stated  that  the  dimensions  for  optimum  fin  geometry  can 
be  expressed  by 
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Figure  4.5  Optimum  Profile  Comparison  Using  Various  Materials 
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where  P^TR),  P2(TR),  and  P3(TR)  are  defined  as  polynomials  in  which  the  taper  ratio 
TR  is  the  independent  variable. 

In  general  the  polynomial  P  is  expressed  in  the  form  of 


n-2 


P(TR)  =  q  TR  "-1  +  C2  77?  n"2  +       +  Cn  TR 


(4.7) 


where  n  is  the  degree  of  the  polynomial. 

The  polynomials  Pl5  P2,  and  P3  are  determined  from  running  the  PROFILE 
AREA  Analysis  Program  using  a  slightly  modified  code.    The  computer  code  listing 
appears  in  Appendix  B.    A  random  number  generator  is  used  to  generate  the  input 
parameters.    The  random  number  generator  produces  input  parameters  in  the  range 
listed  in  Table  12.   Murray-Gardner  (Table  1)  assumption  number  five  (the  fin  width 
is  small  compared  with  its  depth  and  length)  is  checked  in  order  to  insure  that  the 
original  assumptions  are  kept  valid. 
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TABLE  12 


RANGE  OF  THE  RANDOM  INPUT  PARAMETERS 


0.01  <  €  <  0.99 

20  <  L  (  cm  )  <  220 

10  <  k  (  watt  /  m  K)  <  210 

100  <  Tb  (  K  )  <  600 

25  <  qb  (  watts  )  <  1025 

0.01  <  Taper  Ratio  TR  <  0.99 

The  program  enters  into  a  loop  calculating  the  heat  flow  qb  using  the  random 
input  parameters.    The  program  remains  in  this  loop  until  the  maximum  heat  flow  qb 
is  found.    Once  found  the  program  writes  the  random  input  parameters  and  the  value 
determined  to  be  the  maximum  heat  flow  qb  to  a  disk  file.    Another  set  of  random 
input  parameters  is  generated  and  the  maximum  heat  flow  qb  once  again  is  calculated. 
This  sequence  continues  until  100  optimum  data  sets  are  written  to  a  disk  file. 
Three  data  sets  using  the  random  input  parameters  are  produced,  one  each  for  the 
longitudinal  fins  of  rectangular,  trapezoidal,  and  triangular  profiles. 

The  last  step  in  determining  the  polynomials  P,,  P2,  and  P3  is  to  fit  the  collected 
data  stored  in  the  disk  file  to  a  "best"  polynomial  fit  having  the  smallest  degree  d 
using  the  least  square  principle.    Equations  (4.4)  through  (4.6)  are  rewritten  as 
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Lba  €  7\4 
P(TR)X  =  °- *- 


(4.8) 


P(TR),  =  5.    k  a  e  7Y 


/  r  \2 


«*. 


(4.9) 


P(77?)3  =  X/)o  /:  a2  e2  Tfc9 


(4.10) 


Equations  (4.8)  through  (4.10)  are  placed  into  a  MATLAB  "m"  file.    Each  of  the 
collected  data  files  are  read  into  MATLAB  and  were  evaluated  using  the  MATLAB 
polynomial  fitting  function  polyfit  with  the  taper  ratio  TR  as  the  independent  variable 
using  varying  degrees  d  in  order  to  find  the  polynomials  Pl9  P2,  and  P3  which  "best" 
approximate  the  previously  collected  data  from  the  disk  file.   The  results  appear  in 
Table  13.   The  listing  of  the  MATLAB  evaluation  m  files  appear  in  Appendix  C. 
The  standard  error  c  is  defined  as 


£.   = 


1     1/2 


E  fr  -  >',) 2 

N  -  1 


where  %  is  the  computed  value  of  y;  and  there  are  N  values. 


(4.11) 
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TABLE  13 


COEFFICIENTS  OF  THE  OPTIMUM  POLYNOMIALS 


Profile 

Degree 
d 

C, 

C* 

c3 

Standard 
Error 

PI  RECT 
PI  TRAP 
PI  TRI 

2 
2 
2 

0.8675 

-0.0290 

0.0000 

0.0000 
0.8138 
0.9040 

0.0000 
0.0000 
0.0000 

3.54e-3 
2.12e-3 
2.55e-3 

P2  RECT 
P2TRAP 
P2TRI 

2 
3 
2 

1.8648 
0.3881 
0.0000 

0.0000 

-1.0324 

2.4852 

0.0000 
2.7522 
0.0000 

1.63e-2 
1.22e-2 
1.40e-2 

P3  RECT 
P3  TRAP 
P3TRI 

2 
2 
2 

1.6178 
0.5147 
0.0000 

0.0000 
1.1382 
1.1234 

0.0000 
0.0000 
0.0000 

2.06e-2 
1.24e-2 
9.46e-3 

The  equation  for  the  optimum  profile  area  A^  for  either  the  rectangular, 
triangular,  or  trapezoidal  profiles  is  written  as 


A^  =  (  1.1382  +  0.5147  TR  ) 


1 


ko2  e2 


4b 


\3 


L  T 


(4.12) 
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V.    CONCLUSIONS 

This  thesis  has  described  a  new  technique  that  uses  a  parameterization  called  the 
thermal  transmission  matrix  to  represent  a  single  fin  to  evaluate  the  performance  of 
several  fins.   The  algorithm  is  based  on  the  simple  ideal  of  subdividing  the  fin  into 
100  subfins  and  then  the  thermal  transmission  matrix  is  calculated  for  each  subfin. 
Then,  by  a  simple  matrix  multiplication  of  the  individual  subfin  thermal  transmission 
matrices,  the  overall  thermal  transmission  matrix  is  computed.    Not  only  will  this 
produce  the  base  and  tip  conditions  of  the  overall  fin  entity,  but  the  temperature 
distribution  and  heat  flow  at  any  point  between  the  fin  base  and  the  fin  tip  can  be 
easily  determined. 

Several  conclusions  can  be  made  from  the  Chapter  4  Program  Analysis  Results 
and  Discussion.    It  was  observed  that  the  profile  area  Ap,  increases  as  the  cube  of  the 
heat  flow  qb.   In  order  to  double  the  heat  flow  qb,  the  profile  area  Ap  must  be 
increased  eight  times. 

Material  selection  is  very  important.   The  weight  of  the  fin  is  proportional  to 
the  density  and  the  thermal  conductivity  is  inversely  proportional  to  the  profile  area. 
Materials  having  a  thermal  conductivity  greater  than  150  watt  /  m  K  and  a  density 
less  than  3  kg  /  dm3  are  generally  the  best  suited  materials.    For  the  same  material, 
environmental  conditions,  heat  flow,  and  base  temperature,  the  optimum  longitudinal 
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fin  of  triangular  profile  requires  only  69  percent  as  much  material  as  the  optimum 
longitudinal  fin  of  rectangular  profile,  with  the  optimum  longitudinal  fin  of 
trapezoidal  profile  falling  in  between. 

Simple  working  relationships  to  compute  the  optimum  fin  dimensions  are  found. 
Table  16  compares  the  findings  from  this  thesis  with  the  results  from  previous  work 
by  others.    The  differences  between  the  optimum  dimension  polynomial  coefficients 
listed  in  Table  16  are  very  small,  which  verifies  the  accuracy  of  the  thermal 
transmission  matrix  parameterization  technique. 

TABLE  16 
OPTIMUM  DIMENSION  POLYNOMIAL  COEFFICIENTS 


Pi 


P2 


P3 


RECTANGULAR 

.8675 

1.8648 

1.6178 

This  Thesis 

.8846 

1.8473 

1.6341 

Chung  and  Nguyen" 

.8842 

1.8579 

1.6056 

Liu12 

.8846 

1.8473 

1.6341 

Bartas  and  Sellers13 

TRIANGULAR 

.9040 

2.4852 

1.1234 

This  Thesis 

.9485 

2.3104 

1.0957 

Chung  and  Nguyen 

.9546 

2.3027 

1.0991 

Nilson  and  Curry14 

.9091 

2.4255 

1.1025 

Kern  and  Kraus15 

This  new  approach  to  fin  synthesis  and  analysis  permits  an  exact  solution  to  any 
degree  of  accuracy  that  the  user  wishes  to  specify  and  produces  the  result  in  a 
computational  efficient  manner. 
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APPENDIX  A  COMPUTER  CODE  LISTING 


/*  *************************************************************************** 

Main  Program 

****************************************************************************   *  ; 

#include  <process.h> 
#include  <conio.h> 
tfinclude  <stdlib.h> 
#include  <stdio.h> 


#define  Yes 

1 

#define  No 

0 

char 

*S1[] 

=  {  ' 

char 

*S2[] 

=  {  ' 

char 

*S3[] 

=  (  ' 

char 

*S4[] 

=  C  ' 

char 

*S5[] 

=  <:  • 

char 

*S6[] 

=  c  • 

char 

*S7[] 

=  <:  • 

char 

*S8[] 

=  c  • 

char 

*S9[] 

=  <  ' 

watt  /  deg  C",  "  BTU  /  hr  deg  F"  >; 

m  *  2",  "  cm  A  2",  "  ft  '  2",  "  in  *  2",  "  mm  *  2"  }; 
watts",  "  kw",  "  BTU  /  sec",  "  BTU  /  hr"  }; 
deg  F  hr  /  BTU"  }; 
,  "  in",  "  mm"  }; 
Ibm",  "  oz"  }; 
deg  K",  "  deg  C",  "  deg  F",  "  deg  R"  >; 
watt  /  m  deg  K",  "  BTU  /  ft  hr  deg  F"  >; 
No",  "  Yes"); 


deg  C  /  watt",  ' 
m",  "  cm",  "  ft' 
kg",  "  grams",  ' 


char   *0[]   =  {  "<1>  watt  /  deg  C  <2>  BTU  /  hr  deg  F", 

"<1>  m  *  2  <2>  cm  *  2  <3>  ft  *  2  <4>  in  "  2  <5>  mm  *  2", 

"<1>  watts  <2>  kw  <3>  BTU  /  sec  <4>  BTU  /  hr", 

"<1>  deg  C  /  watt  <2>  deg  F  hr  /  BTU", 

"<1>  m  <2>  cm  <3>  ft  <4>  in  <5>  mm", 

"<1>  kg  <2>  grams  <3>  Ibm  <4>  oz", 

"<1>  deg  K  <2>  deg  C  <3>  deg  F  <4>  deg  R", 

"<1>  watt  /  m  deg  K  <2>  BTU  /  ft  hr  deg  F", 

"<1>  No  <2>  Yes"  >; 

void  print_main_menu( ),    print_menu_one( ),    print_menu_two( ); 
void  setupO,    overlaysO,    update_uni  ts(  ); 

float   error=0.1; 

int    T=0,    L=0,    M=0,    S=0,    H=0,    d=4,    K=0,    A=1,    Z=0,    Y=0; 

char  main_key=0; 

int  main()( 

FILE   *data_in; 

if    (      (data_in  =   fopen("uni ts.ext","r"))      !=   NULL)( 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f",&T,&L,&M,&S,&H,&d,&K,&A,&Z,&Y,&error); 
fclose(data_in);    } 

pr  i  nt_ma  i  n_menu( ) ; 

do< 

V  if    (main_key  ==    '1')   setupO; 

if    (main_key  >=    '2'    &&  main_key   <=    '4')   overlaysO;    > 

while   ((   main_key  =   getch()    )    !=    '5'); 
window(1,1,80,25); 
clrscrO; 
_setcursortype(_NORMALCURSOR ) ; 
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return  0;} 

void  over lays(void){ 

char  key=0; 

int  Quit  =  No,  result; 

print_menu_one( ); 

do 

switch  (keyH 

case  'A':  clrscr();  Quit  =  Yes; 

if  (mainjcey  ==  '2')  result  =  spawnl(P_WAIT,  "RECT.EXE",  NULL); 
if  (main_key  ==  '3')  result  =  spawnl(P_WAlT,  "TRAP. EXE",  NULL); 
if  (main_key  ==  '4')  result  =  spawnl(P_WAIT,  "TRI .EXE",  NULL); 
break; 

case  'B':  clrscr();  Quit  =  Yes; 

if  (main_key  ==  '2')  result  =  spawnl (P_WAIT,  "ARECT.EXE",  NULL); 
if  (main_key  ==  '3')   result  =  spawnl (P_WAIT,  "ATRAP.EXE",  NULL); 
if  (main_key  ==  '4')  result  =  spawnl (P_WA1T,  "ATRI.EXE",  NULL); 
break; 

case  'C:  clrscr();  Quit  =  Yes; 

if  (main_key  ==  '2')  result  =  spawnl (PJJAIT,  "0RECT.EXE",  NULL); 

if  (main_key  ==  '3')   result  =  spawnl (P_UAIT,  "OTRAP.EXE",  NULL); 

if  (main_key  ==  '4')  result  =  spawnl (P_WAIT,  "OTRI.EXE",  NULL);} 

while  ((  key  =  toupper(getch( ) )  )  !=  'D'  &&  Quit  ==  No); 

pr  i  nt_ma  i  n_menu( )  ; 
return;} 

void  setup(void){ 
char  key=0; 

print_menu_two( ); 

do 

switch  (keyK 

case  'A':    gotoxy(2, 17);  clreolO;  cprintf ("%s",0[0]  ); 

while  (key  !=  '1'  &&  key  !=  '2')  key  =  getchO; 
Y  =  key  -  ' 1 ' ;print_menu_two( );  break; 

case  'B':    gotoxyO ,  17);  clreolO;  cprintf ("%s",0[1] ); 

while  (key  !=  '1'  &&  key  !=  '2'&S  key  !=  '3'&&  key  !=  'A'&&  key  !=  '5') 

key  =  getch( ); 

A  =  key  -  ' 1 ' ;print_menu_two( );  break; 

case  'C:    gotoxyO,  17);  clreolO;  cprintf  ("%s",0[2]  ); 

while  (key  !=  '1'  &&  key  !=  '2'&&  key  !=  '3'&&  key  !=  '4')  key  =  getchO; 
H  =  key  -  ' 1 ' ;print_menu_two();  break; 

case  'D':    gotoxy(2, 17);  clreolO;  cprintf  ("%s",0[3]  ); 

while  (key  !=  '1'  &&  key  !=  '2')  key  =  getchO; 
Z  =  key  -  ' 1 ' ;print_menu_two();  break; 

case  'E':    gotoxy(2,17);  clreolO;  cprintf  ("%s",0[4]  ); 

while  (key  !=  '1'  &&  key  •=  '2'&&  key  !=  '3'&&  key  !=  '4'&&  key  !=  '5') 

key  =  getchO; 

L  =  key  -  ' 1 ' ;print_menu_two();  break; 

case  'F1:    gotoxy(3, 17);  clreolO;  cprintf ("%s",0[5] ); 

while  (key  !=  '1'  &&  key  !=  '2'&&  key  !=  '3'&&  key  !=  '4')  key  =  getchO; 
M  =  key  -  ■ 1 ' ;print_menu_two();  break; 

case  'G1:    gotoxyO  ,  17);  clreolO;  cprintf  ("%s",0[6]  ); 

while  (key  !=  '1'  &&  key  !=  '2'&&  key  !=  '3'&&  key  !=  '4')  key  =  getchO; 

74 


T  =  key  -  ' 1 ' ;print_menu_two( );  break; 

case  'H':    gotoxyd ,  17);  clreolO;  cprinff  ("%s",0[7] ); 

while  (key  !=  '1'  &&  key  !=  '2')  key  =  getch(); 
K  =  key  -  ' 1 ' ;print_menu_two( );  break; 

case  'I':    gotoxy(2, 17);  clreolO;  cprintf ("%s",0[4] ); 

while  (key  !=  '1'  &&  key  !=  '2'&&  key  !=  '3'&&  key 

key  =  getch( ); 

d  =  key  -  ' 1 ' ;print_menu_two();  break; 

case  'J':    gotoxy(12, 17);  clreolO;  cprintf ("%s",0[8] ); 

while  (key  !=  '1'  &&  key  !=  '2')  key  =  getchO; 
S  =  key  -  ' 1 ' ;print_menu_two();  break; 

case  'K':    gotoxy(22, 12);  clreolO;  scanf ("%f",&error); 
print_menu_two( );> 


•4'&&  key   !=    '5') 


while  ((   key  =  toupper(getchO)    )    !=    'L'); 

update_uni  ts(); 
print_main_menu( ); 
return;) 

void  update_uni ts(void)( 

FILE  *unit; 

unit  =  f openC'uni ts.ext","w"); 

fprintf(unit,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",T,L,M,S(H,d,K,A,Z( Y.error); 

fclose(uni t); 

return;) 

void  print_main_menu( void)( 

_setcursortype(_NOCURSOR ) ; 

wi ndow( 1,1, 80,25 )  ; 

clrscrO; 

gotoxy(27,3);    cprintf ("EXTENDED   SURFACE   ANALYSIS"); 

window(20,6,80,25); 


cprintf ("<1> 
cprintf("<2> 
cprintf ("<3> 
cprintf ("<4> 
cprintf ("<5> 
return;) 


Setup  /  Initial ization\r\n"); 

Fin  of  Rectangular  Prof i le\r\n"); 
Fin  of  Trapezoidal  Prof i le\r\n"); 
Fin  of  Triangular  Prof i le\r\n"); 
Return  to  DOS"); 


void  print_menu_one(voidK 


clrscr( 
gotoxy( 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
return;) 


D; 


■ \r\n"); 

<A>  Predict  the  Performance\r\n"); 
<B>  Analyze  the  Profile  Area\r\n"); 
<C>  Optimize  the  Dimensions\r\n"); 
<D>  Return  to  Main  Menu\r\n"); 

"); 


void  print_menu_two(voidK 


clrscr( ; 

; 

gotoxyC 

,D; 

cprintf  ( 

cprintf  ( 

"<A> 

Admittance 

:  %s\r\n", 

SHY]) 

cprintf  ( 

"<B> 

Area 

:  %s\r\n", 

S2[A]) 

cprintf  ( 

»<c> 

Heat 

:  %s\r\n", 

S3[H]) 

cprintf  ( 

"<D> 

Impedance 

:  %s\r\n", 

S4[Z]) 

cprintf  ( 

»<E> 

Depth 

:  %s\r\n", 

S5[L]) 

\r\n"); 
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cprintf ( 

•<F> 

Mass        : 

%s\r\n". 

S6[M]) 

cprintf ( 

'<G> 

Temperature   : 

%s\r\n", 

S7[T]) 

cprintf ( 

'<H> 

Thermal  K    : 

%s\r\n", 

S8[K]) 

cprintf ( 

'<I> 

Width        : 

%s\r\n", 

S5[d]) 

cprintf ( 

'<J> 

Save  Data    : 

%s\r\n", 

S9[S]) 

cprintf ( 

■<K> 

Percent  Error: 

%f\r\n" 

,  error); 

cprintf ( 

'<L> 

Return  to  Main 

Menu\r\n 

'); 

cprintf ( 

return;} 

"); 
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/*  *************************************************************************** 

PERFORMANCE  Analysis  Program 
Longitudinal  Fin  of  Rectangular  Profile 
****************************************************************************  */ 

^include  <stdio.h> 
#include  <stdlib.h> 
^include  <math.h> 
#include  <conio.h> 
# include  "units.h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  Yb(),  A,  B,  C,  D; 

double  length,  width,  depth; 

double  inc,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data[N+1]  [6] ; 

double  qb=0,A1,  B1,  C1,  D1,  Ta,  Tb,  Ttip,  Tave,  temp,  qa,  segmentjdif f ; 

float  error,  V,  epsilon,  k,  Tbase,  density,  mass,  K2,  eta,  PN; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

void  input(),  computeO,  resultsO; 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main()( 

FILE  *data_in,  *out; 

data_in  =  fopenC'uni ts.ext","r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d   %d  %d  %d   %f ",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

inputO; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  .0001  +  Tave; 

inc  =  .5  *  Ttip;  i  =  1; 

while  (fabs(overall_diff )  >  (Tbase  *  error)   |]  i  <=  N  ){ 

overflow  =  No;      i  =  1;       Ta  =  Ttip;  Tb  =  Ta; 
qa  =  sigma  *  width  *  epsilon  *  pow(Ta,  A)  *  depth; 

while  (i  <=  N  &&  overflow  ==  NoK 

segment_diff  =  10; 

while  (fabs(segment_dif f )  >  ((Tb  *  error)  /  N)  &&  overflow  ==  NoK 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 
temp  =  A  *  Ta  +  B  *  qa; 
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qb       =  C  *   Ta  +  D  *  qa; 
segment_dif f   =   (    temp  -    Tb)   /   2; 
Tb  =   temp  +   segment_dif f ; 

if   (Tb  >  Tbase  +50    ]  ]    Tb  <  OX 

overflow       =  Yes;    possible_overf low  =  Yes; 

if   (Ttip  <  Tbvjnin  )   Tbvjnin  =  Ttip; 

inc  =  (Tbvjnin  -   Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -    inc); 

if   (Ttip  <  Tave)  Ttip  =  Tave;}} 

if   (overflow  ==  NoX 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  A; 
data[i][3]  =  B;  data[i][4]  =  C;  data[i][5]  =  D; 
Ta  =   Tb;    qa  =  qb;    ++i;>> 

if    ( number  jjf_tries  >   29X 

clrscrO;   cprintf ("Fatal   Error:     Unable  to  Compute"); 
getch();    exit(0);) 

if    (overflow  ==   No){ 

if    (Tmax   <   Ttip) 

if    (    Tb  <   Tbase   )   Tmax  =   Ttip; 

overall_diff   =   (Tb  -   Tbase)   /  2; 

if   (   possiblejoverf low  ==  Yes){ 

inc  =   inc  /  2; 

if  (overal l_dif f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal I _d iff; 

numbe  r_of _t r i  es++ ; 

if  (Ttip  <  Tave)  Ttip  =  Tave; 

gotoxy(27,24); 

cprintfC'N   =  %d  Ta   =  %A.1f   Tb  =  %4.1f        ",number_of_tries,w_T(T,Ttip),w_T(T,Tb));» 

for   (i    =  N    -    1;    i    >=   1;    i--)( 

A1  =  A  *  data[i][2]  +  B  *  data[i][4];  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A   =  A1;    B   =   B1;    C=C1;    D=D1;> 

if    (S   ==   1){ 

out  =  fopen(f i le_name,    "w"); 
for    (i    =    1;    i    <=   N;    i++) 

fprintf(out,    "%e,     %e,     %e\n",   w_L(d,    i   *    length  /  N),    w_T(T,data [i]  [1] ),   w_H(H, 
data[i]  [0])    ); 
fclose(out);) 

Ta  =  Ttip; 

mass  =  width  *  length  *  depth  *  density  *  1000; 

PN  =  length  *  sqrt(  sigma  *  epsilon  *  pow(Tbase,3)  /  (  k  *  width)); 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

results( ); 

exit(0); 

return;} 

/ *  *********************************************************************** 

OUTPUT  PROCEDURE 
**************************************************************************  *  / 

void  results(void){ 

clrscrO; 

cprintf("\a  THE   LONGITUDINAL    FIN  OF   RECTANGULAR   PROFILE"); 

cpr  i  ntf ( "\r\n         ===================================================r==============\r\n") ; 
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cpnnt 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
getch( 
return 


Density  =  %4.2f  kg  /  dnT3\r\n",densi ty); 
Fin  Depth  -   ");  P_L(L,  depth); 
Total  Mass  =  ");  p_M(M,  mass); 
Fin  Length  =  ");  p_L(L,  length); 
Fin  Width  =  ");  p_L(d,  width); 
Emissivity  =  %4.2f\r\n",  epsilon); 
Profile  Area  =  ");  p_A(X,  widthMength); 
Profile  Number  -   %4.3f \r\n",PN); 
Base  Heat  Flow  =  ");  p_H(H,  qb); 
Fin  Efficiency  =  %4.3f \r\n",eta); 
Base  Temperature  =  ");  p_T(T,Tb); 
Edge  Temperature  =  ");  p_T(T,Ta); 
Thermal  Conductivity  =  ");  p_K(K,  k); 
Environmental  Parameter  =  %4.2f\r\n",  V); 

Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %4.4f  B  =  ",  A);  p_Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  %A.4f",  D); 


=\r\n"); 


/ 


*  *********************************************************************** 


INPUT  PROCEDURE 


it***********************************************.**************************        *  , 


void  input(void)( 


clrscr( 

);  set 

cprintf 

"\r\n 

cprintf 

"\r\n 

cprintf 

" 

cprintf 

ii 

cprintf 

ii 

cprintf 

H 

cprintf 

ii 

cprintf 

H 

cprintf 

n 

cprintf 

n 

; _set cursor type(_NORMAL CURSOR); 

THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROF  I LE\n\r"); 

Fin  Depth  ");  depth  =  i_L(L); 
Fin  Length  ");  length  -  i_L(L); 
Fin  Width   ");  width  =  i_L(d); 

Emissivity  =  ");  scanf("%f",  Sepsilon); 
Density  (  kg  /  dm"3  )  =  ");  scanf("%f",  &density); 
Base  Temperature  ");  Tbase  =  i_T(T); 
Thermal  Conductivity  ");  k  =  i_K(K); 

Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 


=\r\n"); 


if  (s  ==  D  <: 

cprintf (" 

scanf("%s",  file_name);) 


Output  Filename  is  =  "); 


cprintf ("\r\n=================================================================="); 

if  (  epsilon  >  1.0  ||  epsilon  <=  OX 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getch();  exit(0);> 

if  (  V  <  0  |  ]  V  >  0.8X 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);} 

_setcursortype(_NOCURSOR ) ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  *  / 

double  hr(voidX 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));) 
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double  m(voidX 

return(  sqrt(  (  2.0  *  hr())  /  (  k  *  width)  )  );> 

double  Yb(void){ 

return(  sqrt(  2.0  *  hr()  *  k  *  width)  *  depth  );> 

double  A_(void){ 

return(  cosh(  m( )  *  length  /  N  )   );> 

double  B_(voidX 

return  (  sinh(  m( )  *  length  /  N  )  /  Yb( )  );} 

double  C_(void){ 

return  (  sinh(  m( )  *  length  /  N  )  *  Yb()  );} 

double  D_(void){ 

return(  cosh(  m( )  *  length  /  N  )   );> 
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/*  *************************************************************************** 

PERFORMANCE  Analysis  Program 

Longitudinal  Fin  of  Trapezoidal  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
#include  <stdl ib.h> 
^include  <math.h> 
^include  <conio.h> 
#include  "bessel.h" 
#include  "units. h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  u(double  x),  A,  B,  C,  D; 

double  length,  depth,  edge_width,  base_width,  fin_taper; 

double  b1 ,  db,  a,  b,  ua,  ub,  Ta,  Tb,  a_width,  b_width; 

double  qb=0,  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data[N+1][6] 

double  A1,  B1,  CI,  D1,  taper_ratio,  Ttip,  temp,  qa,  segment_di f f ; 

float  V,  epsilon,  k,  error,  Tbase,  density,  mass,  K2,  Ap,  eta,  PN; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

void  inputO,  resultsO; 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  * / 

int  main(K 

FILE  *data_in,  *out; 

data_in  =  f openC'uni ts.ext", "r"); 

fscanf(data_in,"%d  %d  %d   %d  %d  %d  %d  %d  %d  %d  %f",&T,&L,&M,SS,&H,&d,SK,&X,SZ,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

input ( ); 

taper_ratio  =  edge_width  /  base_width; 
b1  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2*b1  ); 
db  =  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overal l_diff )  >  (Tbase  *  error)   ]J  i  <=  N  K 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  edge_width  *  epsilon  *  pow(Ta,  A)  *  length; 

a  =  b1  -  length; 
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while  (i  <=  H   &&  overflow  ==  NoK 

segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  bl  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No)( 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 

temp  =  A  *  Ta  +  B  *  qa; 

qb   =  C  *  Ta  +  D  *  qa; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_dif f ; 

if   (Tb  >  Tbase  +  50    j  j    Tb  <  OK 

overflow       =   Yes;    possible_overf low  =  Yes; 

if    (Ttip  <   Tbvjnin   )   Tbvjnin  =  Ttip; 

inc   =    (Tbvjnin   -   Tmax)   /   2; 

Ttip  =   fabs(Tbv_min   -    inc); 

if    (Ttip  <   Tave)   Ttip  =  Tave;)) 

if    (overflow  ==  NoH 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i]  [2]  =  A; 
data[i][3]  =  B;  data[i][4]  =  C;  data[i][5]  -  D; 
Ta   =   Tb;    qa   =  qb;    a   =   a  +  db;    ++i;}) 

if    (number_of_tries   >   29){ 

clrscr();    cprintf ("Fatal    Error:      Unable   to  Compute"); 
getchO;    exit(0);} 

if    (overflow  --   No)( 

if    (Tmax   <   Ttip) 

if    (Tb   <   Tbase)    Tmax   =   Ttip; 

overall_diff   =    (Tb   -    Tbase)   /   2; 

if   (   possible_overf low  ==  YesK 

inc   =    inc   /   2; 

if  (overal l_dif f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal Ijdiff; 
number_of _t  r i  es++  ; 
gotoxy(27,24); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f    ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

for    (i    =   N    -    1;    i    >=   1;    i--)( 

A1  -  A  *  data[i][2]  +  B  *  data[i][4];  B1  =  A  *  data [i]  [3]  +  B  *  data[i][5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A  =  A1;    B   =  B1;    C=C1;    D=D1;} 

if    (S  ^    1X 

out   =   fopen(f i le_name,    "w"); 
for    (i    =   1;    i    <=   N;    i++) 

fprintf(out,    "7.e,      %e,      %e\n",    w_L(d,    i    *    length   /   N),    w_T(T,data[i] [1] ),    w_H(H, 
datali]  [0])    ); 
fclose(out);} 

Ta  =  Ttip; 

Ap  =  .5  *  length  *  (  edge_width  +  base_width  ); 

mass  =  Ap  *  length  *  density  *  1000; 

PN  =  2  *  sigma  *  epsilon  *  pow( length,  2)  *  pow(Tbase,3)  /  (k*  base_width); 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  length  *  pow(Tbase,A)); 

results( ); 

exit(0); 

return;} 
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/* 


*.•»«******»*»»******»»*»*«*».«**«*»..»«.*«..».*.».«»*,,».,.,.......... 


OUTPUT  PROCEDURE 


************************************************************************** 


void  results(void){ 


clrscr 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cpr  i  nt 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
getch( 
return 


\r\n 


THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE"); 


Density  =  %4, 

Fin  Depth  =  "), 

Total  Mass  =  "), 

Fin  Length  =  "), 

Emissivity  =  °/A 

Prof  i le  Area  =  "); 

Prof i le  Number  =  %4 

Base  Heat  Flow  =  "), 

Fin  Efficiency  =  %A, 

Fin  Edge  Uidth  =  "); 

Fin  Base  Width  =  "); 

Base  Temperature  =  "); 

Edge  Temperature  =  "), 

Thermal  Conductivity  =  "), 

Environmental  Parameter  =  %A 


2f  kg  /  dm*3\r\n", density); 

p_L(L,  depth); 

p_M(M,  mass); 

p_L(L,  length); 
2f\r\n",  epsilon); 

p_A(X,Ap); 
3f\r\n",PN); 

p_H(H,  qb); 
3f\r\n",eta); 

p_L(d,  edge_width); 

p_L(d,  base  width); 

p_T(T,Tb); 

p_T(T,Ta)| 

p_K(K,  k); 
2f\r\n",  V); 


===========\r\n"); 


Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %A.4f  B  =  ",  A);  p  Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  yj*.UV,    D); 


=\r\n"); 


f* 


****************************************************************** 


INPUT  PROCEDURE 


**************************************************************************  *, 


void  input(void){ 


clrscr( 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 


;  _setcursortype(_NORMALCURSOR); 

"\r\n  THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROF  I LE\r\n"); 

»\r\n    ============================================================r==: 

"  Fin  Depth  ");  depth  =  i_L(L); 

"  Fin  Length  ");  length  =  i_L(L); 

"  Fin  Edge  Width  ");  edge_width  =  i_L(d); 

"  Fin  Base  Width   ");  base_width  =  i_L(d); 

"  Emissivity  =  ");  scanf("%f",  &epsilon); 

"  Density  (  kg  /  dm*3  )  =  ");  scanf("%f",  &density); 

"  Base  Temperature  ");  Tbase  =  i_T(T); 

"    Thermal  Conductivity  ");  k  =  i_K(K); 

"        Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 


=  \r\n"); 


if  (S  ==  1)  { 

cprintf  (" 

scanf("%s",  file_name);> 


Output  Filename  is  =  "); 


cprintf ("\r\n 


="); 


if  (  edge_width  >=  base_width){ 

cprintf("\r\n\r\nFATAL  ERROR: 
getch();  exit(0);) 


Fin  Edge  Width  >=  Fin  Base  Width"); 


if  (  epsilon  >  1.0  ]]  epsilon  <=  0)( 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getch();  exit(0);) 
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if  (  V  <  0  !  !  V  >  0.8K 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getch();  exit(0);} 

_setcursortype(  JJOCURSOR ) ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void)( 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));> 

double  m(voidX 

return(  hr()  /  (k  *  sin(f in_taper)));} 

double  u(double  xX 

return(  2  *  sqrt(x  *  m( )  )  );) 

double  A_(voidX 

return(  ua  *  ( 1 1 (ua)*K0(ub)  +  lO(ub)*Kl (ua) )   );> 

double  B_(void){ 

return  (  ub*ub*( 10(ub)*K0(ua)- I0(ua)*K0(ub) )  /  (  2*k*b_width*depth*m()));} 

double  C_(void){ 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*K1 (ua)  -  1 1 (ua)*K1(ub))/(ua)  );> 

double  D_(void)t 

return(  ub*(  I0(ua)*K1 (ub)+I 1 (ub)*K0(ua))  );) 
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/*  *************************************************************************** 

PERFORMANCE  Analysis  Program 
Longitudinal  Fin  of  Triangular  Profile 

****************************************************************************  » / 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <conio.h> 
#include  "bessel.h" 
^include  "units. h" 

#define  Yes  1 

#define  No  0 

#define  N  99 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  A_(),  B_(),  C_(),  D_(),  m( ) ,  hr(),  u(double  x),  A,  B,  C,  D; 

double  length,  depth,  edge_width,  base_width,  fin_taper; 

double  b1 ,  db,  a,  ua,  ub,  Ta,  Tb,  a_width,  b_width; 

double  qb=0,  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data  [N+1]  [6] 

double  A1,  B1,  C1,  D1,  Ttip,  temp,  qa,  segmentjdi f f ; 

float  V,  epsilon,  k,  error,  Tbase,  density,  mass,  K2,  Ap,  PN,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

void  input(),  resultsO; 


*  *********************************************************************** 


MAIN  PROGRAM 
**************************************************************************  */ 

int  main()( 

FILE   *data_in,  *out; 

data_in  =  f open("uni ts.ext", "r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

input ( ); 

db  =  length  /  (  N  +  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*b1  ); 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  1=1; 

while  (fabs(overal I _d iff )  >  (  Tbase  *  error)   j j  i  <=  N  K 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No)( 
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ub  =  u(db); 

temp  =  Ta  /  IO(ub); 

C    =  2  *  k  *  edge_width  *  depth  *  pow(m(),2)  *  M(ub)  /  (ub  *  I0(ub)); 

qb   =  C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment  jdi ff; 

if    (Tb  >   Tbase  +   50    j  j    Tb  <   OX 

overflow       =   Yes;    possible_overf low  =   Yes; 

if    (Ttip  <   Tbvjnin   )   Tbv_min  =   Ttip; 

inc   =   (Tbvjnin   -    Tmax)   /   2; 

Ttip  =   fabs(Tbv_min  -    inc); 

if    (Ttip  <   Tave)   Ttip  =   Tave;}} 

if   (overflow  ==  NoX 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  1; 
data[i][3]  =  0;  data[i][4]  =  C;  data[i][5]  =  1  ; 
Ta  =  Tb;    qa  =  qb;    ++i;> 

a  =  b1    -    length; 

while  (i  <=  N  &&  overflow  ==  NoK 

segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  b1  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

ua   -   u(a)    ;    ub  =   u(a  +   db); 

A   =   A_();    B   =   B_();    C   =   C_();    D   =   D_(); 

temp  =A*Ta+B*qa; 

qb       =   C  *   Ta  +  D   *  qa; 

segment_diff   =    (    temp   -    Tb)   /   2; 

Tb  =   temp  +   segment  jdi ff; 

if    (Tb  >   Tbase  +  50    ||    Tb  <  OX 

overflow       =   Yes;    possible_overf low  =   Yes; 

if    (Ttip   <   Tbvjnin    )    Tbvjnin  =   Ttip; 

inc   =    (Tbvjnin   -    Tmax)   /   2; 

Ttip  =    fabs( Tbvjnin   -    inc); 

if    (Ttip  <   Tave)   Ttip  =   Tave;)} 

if    (overflow  =■-   NoX 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  A; 
data[i][3]  =  B;  data[i][4]  =  C;  data[i][5]  =  D; 
Ta   =   Tb;    qa  =  qb;    a  =   a  +  db;    ++i;}> 

if   ( number jjf_tries  >  29){ 

clrscrO;   cprintf ("Fatal   Error:     Unable  to  Compute"); 
getch();    exit(0);) 

if    (overflow  ==  NoX 

if    (Tmax   <   Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overalljdiff  =  (Tb  -  Tbase)  /  2; 

if  (  possiblejjverf low  ==  YesX 

inc  =  inc  /  2; 

if  (overal  I  jiif f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;) 
else 

Ttip  =  Ttip  -  overalljdiff; 
number  jjf_t  r  i  es++ ; 
gotoxy(27,24); 

cprintf("N  =  %d  Ta  =  %A.1f  Tb  =  %A.1f   ",numberjjf_tries#w_T(T,Ttip),w_T(T,Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 
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for    (i    =   N    -    1;    i    >=    1;    i--)< 

A1  =  A  *  data[i][2]  +  B  *  data [i]  [4];  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
CI  =  C  *  datatil  [2]  +  D  *  data[i][4];  D1  =  C  *  data [i] [3]  +  D  *  data[i][5]; 
A   =  A1;    B   =  B1;    C=C1;    D=D1;} 


length   =    length   +  db; 

if    (S  ==   1X 

out   =   f open(f i le_name,    "w"); 

for    (i    =    1;    i    <=   N;    i++) 

fprintf(out,    "%e,     %e, 
datati]  [0])); 

fclose(out);} 


%e\n",    w_L(d,    i   *    length  /  N),   w_T(T,data[i] [1] ),    w_H(H( 


Ta   =   Ttip; 

Ap  =    .5   *    length   *   base_width; 

mass   =  Ap  *   depth  *   density  *    1000; 

PN  =  2  *  sigma  *  epsilon  *  pow( length, 2)   *  pow(Tbase,3)   /   (k*base_width); 

eta  =  qb  /   (   2  *  sigma  *  epsilon  *   length  *  depth  *  pow(Tbase,4)); 

result s( ); 

exit(0); 

return;} 


*   *********************************************************************** 


OUTPUT    PROCEDURE 


**************************************************************************    *  , 


void   resul ts(void){ 


clrscr 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
getch( 
return 


'\a 


THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROF  I LE\r\n") ; 


Densi  ty 

Fin  Depth 

Total  Mass 

Fin  Length 

Emissivity 

Prof i le  Area 

Prof i le  Number 

Base  Heat  Flow 

Fin  Efficiency 

Fin  Base  Width 

Base  Temperature 

Edge  Temperature 

Thermal  Conductivity 

Environmental  Parameter 


%4.2f  kg  /  dmA3\r\n", density); 

");  p  L(L,  depth); 

");  p_M(M,  mass); 

");  P_L(L,  length); 

%A.2f\r\n",  epsilon); 

");  p_A(X,Ap); 

%4.3f\r\n",PN); 

");  p  H(H,  qb); 

%4.3f\r\n",eta); 

");  p_L(d,  base_width); 

");  p_T(T,Tb);  " 

");  p_T(T,Ta); 

");  p_K(K,  k); 

%4.2f\r\n",  V); 


=\r\n»); 


Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %4.4f  B  =  ",  A);  p_Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  %A.4f",  D); 


=\r\n»>; 


/*  *********************************************************************** 


INPUT  PROCEDURE 


**************************************************************************  * , 


void  input(void)< 


clrscrO;   _setcursortype(_NORMALCURSOR); 

cprintf("\r\n  THE    LONGITUDINAL    FIN   OF   TRIANGULAR   PROFILE\r\n"); 

cprintf ("\r\n         ==============================================r ======== 

cprintf("  Fin  Depth  ");   depth  =  i_L(L); 

cprintf("  Fin  Length  ");  length  =  i  L(L); 


=\r\n»); 
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cprintf("  Fin  Base  Uidth  ");  base_width  =  i_L(d); 

cprintf  ("  Emissivity  =  ");  scanf("%f",  fiepsilon); 

cprintf ("  Density  (  kg  /  dm*3  )  =  ");  scanf("%f",  &density); 

cprintf("  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintf("  Thermal  Conductivity  ");  k  =  i_K(K); 

cprintf ("  Environmental  Parameter  (0  to  .8)  =  ");  scanf("/£f",  &V); 

if  (S  ==  1)  { 

cprintf("  Output  Filename  is  =  "); 

scanf("%s",  file_name);} 

Cprintf("\r\n      ============r==========r===============r=rrr=================r=r=="); 

if  (  epsilon  >  1.0  |J  epsilon  <=  OK 

cprintf ("\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getch();  exit(0);) 

if  (  V  <  0  j  i  V  >  0.8K 

cprintf ("\r\n\r\nFATAL  ERROR:   Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);) 

edge_width  =  base_width  /  1e2; 
_setcursortype(_NOCURSOR ) ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));} 

double  m(void){ 

return(  hr()  /  (k  *  sin(f in_taper)) );} 

double  u(double  x){ 

return(  2  *  sqrt(x  *  m( )  )  );} 

double  A_(void)( 

return(   ua  *  (I 1(ua)*K0(ub)  +  I0(ub)*Kl (ua) )   );} 

double  B_(void){ 

return  (   ub*ub*( I0(ub)*K0(ua)- l0(ua)*K0(ub))  /  (  2*k*b_width*depth*m( )));} 

double  C_(voidK 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*K1 (ua)  -  1 1 (ua)*Kl (ub) )/(ua)  );} 

double  D_(void){ 

return(  ub*(  I0(ua)*K1 (ub)+I 1 (ub)*K0(ua))  );} 
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/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program 
Longitudinal  Fin  of  Rectangular  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
^include  <stdl ib.h> 
# include  <math.h> 
#include  <conio.h> 
#include  "units. h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  length,  width,  depth,  starting_width,  ending_width,  width_inc; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  computet); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(){ 

FILE  *data_in,  *out; 

data_in  =  fopen("uni ts.ext","r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

input ( ); 

width_inc  -    (ending_width  -  start ing_width  )  /  100; 
width  =  starting_width; 

while  (width  <  ending_width)( 

computet ); 

while  (qb  ==  0)( 

width  =  1.001  *  width; 
computeC );> 

gotoxyC 19,23) ; 

cprintfC'width  =  %4.1f        length  =  %4.1f   heat  =  %4.1f",    w_L(d, width),    w_L(L, length), 
w_H(H,    qb)); 

if  (S  ==  1){ 

out  =  f open(f i le_name,  "a"); 

fprintf(out,  "%e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e\n",  w_L(d, width), 

w_T(T,Tb),  w_T(T,Ta),  eta,  w_H(H,  qb),  w_M(M,mass),  A,  w_Z(Z,B),  w_Y(Y,C),  D); 
fclose(out );) 

width   =   width   +   width_inc;> 

gotoxy(21 ,24);clreol( );cprintf ("\aCompleted  \r\n"); 

clreolO; 

getchO; 
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exit(O); 
return  0;} 


/*  *********************************************************************** 


COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(voidX 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  Yb(); 

double  inc,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data[N+1] [6] ; 

double  A1,  B1,  CI,  D1,  Ttip,  temp,  qa,  Tave,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

length  =  Ap  /  width; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  -    .0001  +  Tave; 

inc  =  .5  *  Ttip;  i  =  1; 

while  (f abs(overal l_di f f )  >  (  Tbase  *  error)   j j  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  width  *  epsilon  *  pow(Ta,  4)  *  depth; 

while  (i  <=  N  &&  overflow  ==  NoX 

segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoX 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 
temp  =  A  *  Ta  +  B  *  qa; 
qb   =  C  *  Ta  +  D  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +  50  j  j  Tb  <  OX 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbvjnin  )  Tbv_min  =  Ttip; 

inc  =  (Tbvjnin  -  Tmax)  /  2; 

Ttip  -   fabs( Tbvjnin  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

if    (overflow  ==   No){ 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  A; 
data[i][3]  =  B;  data[i] [4]  =  C;  datati] [5]  =  D; 
Ta   =   Tb;    qa   =  qb;    ++i;)> 

if    ( number jjf_tries   >   29 )( 
qb=0; return;) 

if   (overflow  ==  No){ 

if    (Tmax   <   Ttip) 

if    (Tb  <   Tbase   )   Tmax  =   Ttip; 

overalljdiff  =   (Tb  -   Tbase)  /  2; 

if   (   possiblejjverf low  ==  YesX 
inc   =    inc   /   2; 
if    (overalljdiff   >   0    )      Ttip  =   Ttip   -    inc; 
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else     Ttip=Ttip+inc;} 
else 

Ttip  =   Ttip   -    overal l_dif f  ; 
number_of_tries++; 
gotoxy(27,24); 

cprintfCN   =  %d  Ta  =  %4.1f   Tb  =  %4.1f        ",number_of_tries,w_T(T,Ttip),w_T(T(Tb)); 
if   (Ttip  <  Tave)   Ttip  =  Tave;)) 

for   (i    =   N    -    1;    i    >=    1;    i--K 

A1  =  A  *  data  [i]  [2]  +  B  *  data  [i]  [4];  B1  =  A  *  data  [i]  [3]  +  B  *  data  [i]  [5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i]  [3]  +  D  *  data[i][5]; 
A   =  A1;    B   =   B 1 ;    C=C1;    D=D1;) 

Ta  =  Ttip; 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

mass  =  width  *  length  *  depth  *  density  *  1000; 

return;} 

z*  *********************************************************************** 

INPUT  PROCEDURE 
**************************************************************************  */ 

void  input (voidK 

clrscr( );  _setcursortype(_NORMALCURSOR); 

cprintf("\r\n         THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE"); 

cprintf ("\r\n    =========== =====  =  =  =  =  =  =  ===  =  ===  ===  ===  ======== ======================= \r\n"); 

cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Starting  Fin  Width  ");  start ing_width  =  i_L(d); 

cprintfC  Ending  Fin  Width  ");  ending_width  =  i_L(d); 

cprintfC  Profile  Area  ");  Ap  =  i_A(X); 

cprintfC  Emissivity  =  ");  scanf("%f",  Sepsilon); 

cprintfC  Density  (  kg  /  dm"3  )  =  ");  scanf("%f",  Sdensity); 

cprintfC  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintfC    Thermal  Conductivity  ");  k  =  i_K(K); 

cprintfC         Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 

if    (S   ==   1)    ( 

cprintfC  Output   Filename   is  =  ");    scanf("%s",    file_name);> 

cprintf("\n         ==================================================================\r\n"); 

if  (  epsilon  >  1.0  ||  epsilon  <=  0){ 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getch();  exit(O);) 

if  (  V  <  0  ]  |  V  >  0.8X 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);} 

_setcursortype(_NOCURSOR ) ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));} 

double  m(void)( 

return(  sqrt(  (  2.0  *  hr())  /  (  k  *  width)  )  );} 

double  Yb(void){ 

91 


return(  sqrt(  2.0  *  hr()  *  k  *  width)  *  depth  );> 

double  A_(voidK 

return(  cosh(  m( )  *  length  /  N  )   );} 

double  B_(void){ 

return  (  sinh(  m()  *  length  /  N  )  /  Yb( )  );} 

double  C_(void)C 

return  (  sinh(  m()  *  length  /  N  )  *  Yb()  );} 

double  D_(void){ 

return(  cosh(  m( )  *  length  /  N  )   );) 
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/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program 
Longitudinal  Fin  of  Trapezoidal  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <conio.h> 
#include  "units. h" 
#include  "bessel .h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  length,  edge_width,  base_width,  a_width,  b_width,  depth,  start ing_width,  ending_width; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a,  width_inc; 

float  V,  K2,  epsilon,  k,  Tbase,  taper_ratio,  density,  mass,   Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  computed); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main( ){ 

FILE  *data_in,  *out; 

data_in  =  fopen("uni ts.ext", "r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L ,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

f close(data_in); 

error  =  error  /  100; 

input( ); 

width_inc  =  (ending_width  -  start ing_width  )  /  1 00; 
base_width  =  start ing_width; 

while  (base_width  <  ending_width)( 

compute( ); 

while  (qb  ==  OK 

base_width  -   1.001  *  base_width; 
compute( );) 

gotoxy( 19,23 ); 

cprintf("width  =  %A.1f   length  =  %4.1f  heat  =  %A.1f",  w_L(d,base_width),  w_L(L, length), 
w_H(H,  qb)); 

if    (S   ==   1){ 

out   =   fopen(f i le_name,    "a"); 

fprintf(out,    "%e,     %e,      %e,      %e,     %e,     %e,     %e,     %e,     %e,     %e\n",    w_L(d,base_width), 
w_T(T,Tb),    w_T(T,Ta),    eta,    w_H(H,    qb),    w_M(M,mass),    A,    w_Z(Z,B),    w_Y(Y,C),   D); 
f close(out );> 

base_width  =  base_width  +  width_inc;> 

gotoxy(21,24);clreol();cprintf ("\aCompleted  \r\n"); 
clreolO; 
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getchO; 
exit(O); 
return  0;} 


/*  *********************************************************************** 


COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(voidX 

double  A_(),  B_(),  C_(),  D_( ) ,  m(),  hr(),  u(double  x); 

double  length_inc,  bl,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000,  data[N+1] [6] ; 

double  A1,  B1,  C1,  D1,  Ttip,  temp,  qa,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

edge_width  =  base_width  *  taper_ratio; 
length  =  2  *  Ap  /  (  base_width  +  edge_width  ); 
b1  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2*b1  ); 
db  =  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase, A); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i  =  1; 

while  (f abs(overal l_di f f )  >  (  Tbase  *  error)  \\    i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  edge_width  *  epsilon  *  pow(Ta,  4)  *  depth; 

a  =  b1  -  length; 

while  (i  <=  N  &&  overflow  ==  No)( 

segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  b1  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 

temp  =A*Ta+B*qa; 

qb   =  C  *  Ta  +  D  *  qa; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  j  J  Tb  <  0)( 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =    fabs(Tbv_min   -    inc); 

if    (Ttip  <   Tave)    Ttip  =   Tave;)} 

if    (overflow  ==   NoK 

data[i][0]  =  qb;  data[i] [1]  =  Tb;  dataCi] [2]  =  A; 
data[i][3]  =  B;  data[i]  [4]  =  C;  data[i]  [5]  =  D; 
Ta  =  Tb;   qa  =  qb;   a  =  a  +  db;   ++i;)) 

if    (number_of_tries   >   29){ 
qb=0;    return;) 
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if  (overflow  ==  No)( 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overalt_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf Low  ==  Yes){ 

inc  =  inc  /  2; 

if    (overal l_dif f   >   0    )      Ttip  =   Ttip   -    inc; 

else     Tt ip  =  Ttip  +   inc;) 
else 

Ttip  =   Ttip   -    overal l_di ff; 
number_of _t  r i  es++ ; 
gotoxy( 27,24 ); 

cprintfCN   =  %d  Ta   =  %4.1f   Tb  =  %4.1f      ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)); 
if   (Ttip  <  Tave)   Ttip  =  Tave;}) 

for   (i    =  N   -    1;    i   >=   1;    i--)( 

A1  =  A  *  data[i][2]  +  B  *  data  [i]  [4] ;  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A   =  A1;    B   =   B1;    C=C1;    D=D1;> 

Ta  =  Ttip; 

mass  =  .5  *  (  edge_width  +  base_width  )  *  length  *  depth  *  density  *  1000; 
eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,  4)); 
return; ) 

/*  *********************************************************************** 

INPUT  PROCEDURE 
**************************************************************************  */ 

void  input (voidH 

c I rscr( ) ;  _set cursor type(_NORMAL CURSOR ) ; 

cprintfC  THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE"); 

cprintf("\r\n    ====================== ==========================================r=\r\n»); 

cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Starting  Base  Fin  Width  ");  start ing_width  =  i_L(d); 

cprintfC  Ending  Base  Fin  Width  ");  ending_width  =  i_L(d); 

cprintf("  Taper  Ratio  =  ");  scanf("%f",  &taper_ratio); 

cprintfC  Profile  Area  ");  Ap  =  i_A(X); 

cprintfC  Emissivity  =  ");  scanf("%f",  Sepsilon); 

cprintfC  Density  (  kg  /  dm'3  )  =  ");  scanf("%f",  &density); 

cprintfC  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintfC    Thermal  Conductivity  ");  k  =  i_K(K); 

cprintfC  Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 

if  (S  ==  1)  { 

cprintfC  Output  Filename  is  =  ");  scanf("%s",  file_name);) 

Cprintf("\n      =========r=r=r=rrrr=rrrrr===r===rr================================\r\n"); 

if  (  taper_ratio  >  .99  jj  taper_ratio  <=  0.0)  { 

cprintf ("\r\n\r\nFATAL  ERROR:   Taper  Ratio  is  not  between  0  and  1"); 
getchO;  exit(O);) 

if  (  start ing_width  >=  ending_width){ 

cprintf C\r\n\r\nFATAL  ERROR:   Starting  Base  Fin  Width  >=  Ending  Base  Fin  Width"); 
getchO;  exit(0);) 

if  (  epsilon  >  1.0  ]]  epsilon  <=  0)( 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getchO;  exit(0);) 

if  (  V  <  0  | J  V  >  0.8)( 

cprintf ("\r\n\r\nFATAL  ERROR:   Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);} 
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_setcursortype(_NOCURSOR ) ; 
return;) 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));> 

double  m(void)( 

return(  hr()  /  (k  *  sin(f in_taper)));> 

double  u(double  x){ 

return(  2  *  sqrt(x  *  m( )  )  );) 

double  A_(void){ 

returnC  ua  *  (I1(ua)*K0(ub)  +  lO(ub)*K1(ua))  );} 

double  B_(void){ 

return  (  ub*ub*( I0(ub)*K0(ua)- I0(ua)*K0(ub) )  /  (  2*k*b_width*depth*m( )));} 

double  C_(void){ 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*KHua)  -  11 (ua)*K1 (ub))/(ua)  );} 

double  D_(void){ 

return(  ub*(  I0(ua)*Kl (ub)+I 1 (ub)*K0(ua) )  );> 
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/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program 

Longitudinal  Fin  of  Triangular  Profile 
***************************************************************************  */ 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <conio.h> 
# include  "units.h" 
#include  "bessel .h" 

tfdefine  Yes  1 

tfdefine  No  0 

#define  N  99 

#define  sigma  5. 66961  e- 8 

char  f i  le_name[20] ; 

double  length,  width_edge,  base_width,  a_width,  b_width,  depth,  start ing_width,  ending_width; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a,  inc_width; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  computet); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(X 

FILE  *data_in,  *out; 

data_in  =  fopen("uni ts.ext","r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

inputO; 

inc_width  =  (ending_width  -  start ing_width  )  /  100; 
base_width  =  starting_width; 

while  (base_width  <  ending_width){ 

compute( ); 

while  (qb  ==  0){ 

base_width  =  1.001  *  base_width; 
compute( );} 

gotoxy(19,23);  cprintf ("width  =  %A.1f   length  =  %4.1f  heat  =  %4.1f",  w_L(d,base_width), 

w_L(L, length),  w_H(H,  qb)); 

if  (S  ==    1){ 
out  =  fopen(f i le_name,  "a"); 

fprintf(out,    "%e,     %e,     %e,      %e,     %e,     %e,     %e,     %e,     %e,     %e\n",   w_L(d,base_width), 
w_T(T,Tb),    w_T(T,Ta),    eta,    w_H(H,    qb),    w_M(M,mass),    A,    w_Z(Z,B),    w_Y(Y,C),   D); 
fclose(out);) 

base_width  =  base_width  +   inc_width;} 

gotoxy(21,2A);clreol( ); cprintf ("\aCompleted  \r\n"); 

clreolO; 

getch(); 
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exit(O); 
return  0;) 


/*  ******** *************************  ******************* ******************* 


COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(voidK 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  u(double  x); 

double  length_inc,  b1 ,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data[N+1] [6] ; 

double  A1,  B1,  CI,  D1,  Ttip,  temp,  qb,  segmentjdif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

length  =  2  *  Ap  /  base_width; 

width  edge  =  base  width  /  100; 

db  =  length  /  <  N~+  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*b1  >; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overal l_di f f )  >  (  Tbase  *  error)  \\    i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No){ 

ub  =  u(db); 

temp  =  Ta  /  I0(ub); 

C    =  2  *  k  *  width_edge  *  depth  *  pow(n(),2)  *  I1(ub)  /  (ub  *  I0(ub)); 

qb   =  C  *  Tb; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if    (Tb  >   Tbase  +   50    j  j    Tb  <   OK 

overflow       =   Yes;    possible_overf low  =   Yes; 

if    (Ttip  <   Tbvjnin    )    Tbv_min   =   Ttip; 

inc   =    (Tbvjnin   -    Tmax)   /   2; 

Ttip  =   fabs(Tbvjnin   -    inc); 

if    (Ttip  <   Tave)    Ttip  =   Tave;}) 

if    (overflow  ==   No)( 

data[i][0]    =  qb;    data[i][1]    =   Tb;    data[i][2]    =   1; 
data[i][3]    =  0;     data[i]  [4]    =  C;     data[i][5]    =   1; 
a  =   b1    -    length; 
Ta  =  Tb;   qa  =  qb;    ++i ;> 

while  (i    <=  N  &&  overflow  ==  No)( 

segment  jiiff  =   10; 

a_width  =  widthj^dge  *  a  /   (   b1    -    length); 

bjwidth   =   widthjfdge  *   (    a  +  db   )   /   (    b1    -    length); 

while   (fabs(segment   diff)    >    (    Tb  *   error)   /   N  &&  overflow  ==  No){ 
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ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 

temp  -   A  *  Ta  +  B  *  qa; 

qb   =  C  *  Ta  +  D  *  qa; 

segment_dif f  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_dif f ; 

if  (Tb  >  Tbase  +  50  j  j  Tb  <  OH 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;}) 

if    (overflow  ==   No){ 

data[i][0]    =  qb;   data[i][1]    =  Tb;   data[i]  [2]    =  A; 

data[i][3]    =   B;      data[i][4]    =  C;      data[i][5]    =  D; 

Ta  =   Tb;    qa   =  qb;    a  =   a  +  db;  ++i;>> 

if   (number_of_tries  >  29 ){ 
qb=0;    return;) 

if   (overflow  ==  No){ 

if   (Tmax  <  Ttip) 

if    (Tb  <   Tbase)   Tmax   =   Ttip; 

overall_diff   =    (Tb   -    Tbase)   /   2; 

if   (   possible_overf low  ==  Yes){ 

inc   =    inc   /   2; 

if  (overal l_dif f  >  0  )   Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal l_di ff; 
number_of _t  r  i  es++  ; 
gotoxy(27,24); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)} 

for   (i    -  N   -    1;    i    >-   1;    i--)( 

A1  -  A  *  data[i][2]  +  B  *  data[i][A];  B1  -  A  *  data[i][3]  +  B  *  data  [i]  [5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A   =   A1;    B   =   B 1 ;    C=C1;    D=D1;) 

length   =    length   +   db; 

Ta   =   Ttip; 

mass  =  .5  *  base_width   *  length  *  depth  *  density  *  1000; 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4) ); 

return;) 


/ 


*  *********************************************************************** 


INPUT  PROCEDURE 


**************************************************************************  * , 


void  input (voidX 


clrscr( 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 


_setcursortype(_NORMALCURSOR ) ; 

THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROF  I LE\r\n"); 


Ar\n"): 


Fin  Depth  ");  depth  =  i_L(L); 
Starting  Base  Fin  Width  ");  start ing_width  =  i_L(d); 
Ending  Base  Fin  Width  ");  ending_width  =  i_L(d); 
Profile  Area  ");  Ap  =  i_A(X); 

Emissivity  =  ");  scanf("%f",  &epsilon); 
Density  (  kg  /  dnT3  )  =  ");  scanf("%f",  &density); 
Base  Temperature  ");  Tbase  =  i_T(T); 
Thermal  Conductivity  ");  k  =  i_K(K); 
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cprintfC        Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 

if  (S  ==  1)  { 

cprintfC  Output  Filename  is  =  ");  scanf("%s",  f i  I  e_name ) ; } 

cprintf ("\n    ==================================================================\r\n"); 

if  (  start ing_width  >=  ending_width){ 

cprintf ("\r\n\r\nFATAL  ERROR:   Starting  Base  Fin  Width  >=  Ending  Base  Fin  Width"); 
getch();  exit(0);) 

if  (  epsilon  >  1.0  ]|  epsilon  <=  OK 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getch();  exit(0);} 

if  (  V  <  0  | |  V  >  0.8)( 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getch();  exit(0);) 

_setcursortype(_N0CURS0R )  ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));} 

double  m(void){ 

return(  hr()  /  (k  *  sin(f in_taper) ));} 

double  u(double  xK 

return(  2  *  sqrt(x  *  m( )  )  );) 

double  A_(void){ 

return(   ua  *  ( 1 1 (ua)*K0(ub)  +  I0(ub)*Kl (ua) )   );) 

double  B_(void){ 

return  (  ub*ub*( 10(ub)*K0(ua)- I0(ua)*K0(ub))  /  (  2*k*b_width*depth*m( )));} 

double  C_(void){ 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*K1(ua)  -  11 (ua)*K1 (ub))/(ua)  );) 

double  D_(void){ 

return(  ub*(  I0(ua)*K1(ub)+I 1(ub)*K0(ua))  );) 

/*  *************************************************************************** 

OPTIMIZE  Analysis  Program 

Longitudinal  Fin  of  Rectangular  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
# include  <conio.h> 
#include  "units. h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

double       length,  width,  depth,  A,  B,  C,  D,  Ta,  Tb,  qb=0; 
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float        Tbase,  heat,  PN,  eta,  V,  K2,  epsilon,  k,  density,  mass,  Ap,  error; 
int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void         inputO,  compute(),  resultsO; 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  *  / 

int  main()( 

FILE  *data_in; 

data_in  =  fopen("uni ts.ext","r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ',&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  1 00; 

input ( ); 

heat  =  heat  /  (  1  -  V); 

Ap  -    (1.6178  *  pow(heat/depth  ,  3))  /  (  k  *  pow(epsi Ion, 2)  *  pow(Tbase,9)  *  pow(sigma,2)); 

width  =  (1.8648*  pow(heat/depth  ,  2))  /  (  k  *  sigma*epsi lon*pow(Tbase,5) ); 

computet ); 
results( ); 
exit(0); 

return  0;> 

/*  *********************************************************************** 

COMPUTE  ROUTINE 
**************************************************************************  */ 

void  computet voidX 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  Yb(); 

double  inc.  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000,  data[N+1] [6] ; 

double  A1,  B1,  C1,  D1,  Ttip,  Tave,  temp,  qa,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

length  =  Ap  /  width; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 

Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  .0001  t  Tave; 

inc  =  .5  *  Ttip;  i  =  1; 

while  (fabs(overal l_dif f )  >  (  Tbase  *  error  >  J j  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  width  *  epsilon  *  pow(Ta,  4)  *  depth; 

while  (i  <=  N  &&  overflow  ==  NoK 

segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

A  =  A _();  B  =  B _();  C  =  C_();  0  =  D_(); 
temp  =  A  *  Ta  +  B  *  qa; 
qb   =  C  *  Ta  +  0  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
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Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +  50  j  J  Tb  <  OX 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2  ; 

Ttip  =   fabs(Tbv_min   -    inc); 

if    (Ttip  <   Tave)   Ttip  =   Tave;)} 

if   (overflow  ==  NoK 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  A; 
data[i][3]  =  B;  data[i][A]  =  C;  data[i][5]  =  D; 
Ta  =   Tb;    qa  =  qb;    ++i;}} 

if    (number_of_tries   >   29  H 

clrscrO;   cprintf ("Fatal   Error:     Unable  to  Compute"); 
getchO;   exit(0);> 

if   (overflow  ==  NoH 

if    (Tmax   <   Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  --   YesK 

inc  =  inc  /  2; 

if  (overal l_di f f  >0)  Ttip  =  Ttip  -  inc ; 
else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal I _di ff; 
number_of _t  r  i  es++ ; 
gotoxy(27,24); 

cprintfC'N  =  %d  Ta  -   %4.1f  Tb  =  %4.1f    ",number_of_tries,w_T(T,Ttip),w_T(T, Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

for   (i    =   N    -    1;    i    >=   1;    i--)< 

A1  =  A  *  data[i][2]  +  B  *  data[i][4];  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
CI  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A   =   A1;    B   =   B1;    C=C1;    D=D1;> 

Ta  =  Ttip; 

eta  =  qb  /  (2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

PN  =  length  *  sqrt(sigma  *  epsilon  *  pow(Tbase,3)  /  (k  *  width)); 

mass  =  width  *  length  *  depth  *  density  *  1000; 

return;} 


*  *********************************************************************** 


OUTPUT  PROCEDURE 


**************************************************************************  */ 


void  resul ts(void)( 


clrscr( ) 

cprintf ( 

"\a 

cprintf ( 

"\r\n 

cprintf ( 

"\r\n 

cprintf ( 

II 

cprintf ( 

11 

cprintf ( 

II 

cprintf ( 

II 

cprintf ( 

II 

cprintf ( 

11 

cprintf ( 

II 

cprintf ( 

II 

cprintf ( 

II 

cprintf ( 

II 

cprintf ( 

II 

OPTIMUM  PARAMETERS"); 
THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE"); 

Density  =  %A.2f  kg  /  dmA3\r\n",densi ty); 
Fin  Depth  =  ");  p_L(L,  depth); 
Total  Mass  =  ");  p_M(M,  mass); 
Fin  Height  =  ");  p_L(L,  length); 
Fin  Width  =  ");  p_L(d,  width); 
Emissivity  =  %4.2f\r\n",  epsilon); 
Profile  Area  =  ");  p_A(X,  Ap); 
Profile  Number  =  %A.3f \r\n",PN); 
Base  Heat  Flow  =  ");  p_H(H,  qb); 
Fin  Efficiency  =  %4.3f \r\n",eta); 
Base  Temperature  =  ");  p_T(T,Tb); 


=\r\n"); 
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cprintf 

cprintf 

cprintf 

cprintf 

cprintf 

cprintf 

cprintf 

cprintf 

getchO; 

return;} 


Edge  Temperature  =  ");  p_T(T,Ta); 
Thermal  Conductivity  -   ");  p_K(K,  k); 
Environmental  Parameter  =  %4.2f\r\n",  V); 


Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %4.4f  B  =  ",  A);  p_Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  U.Af",  D); 


=\r\n"); 


/*  I****-********-*********************************************************** 


INPUT  PROCEDURE 


**************************************************************************  *, 


void  input(void){ 


clrscr( 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 


setcursortype(_NORMALCURSOR ) ; 
'\r\n  THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE\r\n"); 

i\r\n    ============================================================== 

1  Fin  Depth  ");  depth  =  i_L(L); 

1  Heat  Dissipation  ");  heat  =  i_H(H); 

1  Emissivity  =  ");  scanf("%f",  Sepsilon); 

1  Density  (  kg  /  dnT3  )  =  ");  scanf("%f",  Sdensity); 

1  Base  Temperature  ");  Tbase  =  i_T(T); 

1  Thermal  Conductivity  ");  k  =  i_K(K); 

'  Environmental  Parameter  (0  to  .8)  -   ");  scanf("%f",  &V); 

'\r\n    ==================================================================»); 


\r\n"); 


if  (  epsilon  >  1.0  \\    epsilon  <=  0){ 

cprintf ("\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(0);> 


if  (  V  <  0 


V  >  0.8){ 


cprintf ("\r\n\r\nFATAL  ERROR:   Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);} 

_setcursortype(_NOCURSOR ) ; 
return;} 

/*  *********************************************************************** 

GLOBAL  FUNCTIONS 
**************************************************************************  */ 

double  hr(void)( 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));} 

double  m(void){ 

return(  sqrt(  (  2.0  *  hr())  /  (  k  *  width)  )  );} 

double  Yb(void){ 

return(  sqrt(  2.0  *  hr()  *  k  *  width)  *  depth  );} 

double  A_(void)( 

return(  cosh(  m()  *  length  /  N  )  );} 

double  B_(void)( 

return  (  sinh(  m( )  *  length  /  N  )  /  Yb()  );} 

double  C_(void)( 

return  (  sinh(  m( )  *  length  /  N  )  *  Yb()  );} 

double  D  (void)C 
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return(  cosh(  m( )  *  length  /  N  )   );} 
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/*  *************************************************************************** 

OPTIMIZE  Analysis  Program 

Longitudinal  Fin  of  Trapezoidal  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
# include  <stdlib.h> 
#include  <math.h> 
#include  <conio.h> 
#include  "units. h" 
#include  "bessel .h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5. 66961 e- 8 

char  f i le_name[20] ; 

double  length,  edge_width,  a_width,  b_width,  base_width,  depth,  inc_width; 

double  A,  B,  C,  D,  Ta,  a,  Tb,  qb=0,  ua,  ub,  fin_taper,  Constant_3; 

float  K2,  taper_ratio,  V,  eta,  PN,  epsilon,  heat,  k,  Tbase,  density,  mass,  Ap,  error; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  computet),  resultsC); 

FILE  *data_in,  *out; 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(){ 

data_in  =  f openC'uni ts.ext", "r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L,&M,&S,&H,Sd,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

input ( ); 

heat  =  heat  /  (1  -  V); 

Constant_3  =  0.3881*pow( taper_rat i o, 2)  -  1 .0324*taper_rat io  +  2.7522; 

Ap      =  ((0.5147  *  taper_ratio  +  1.1382)  *  pow(heat/depth  ,  3))  /  (  k  *  pow(epsi Ion, 2)  * 

pow(Tbase,9)  *  pow(sigma,2)); 
base_width  =  (Constant_3  *  pow(heat/depth  ,  2))  /  (  k  *  sigma*epsi lon*pow(Tbase,5)); 
edge_width  -   taper_ratio  *  base_width; 

compute ( ); 
resultsO; 
ex  i  t ( 0 )  ; 
return  0;) 

/*  *********************************************************************** 

COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(void){ 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  u(double  x); 

double  length_inc,  b1 ,  db; 

double  inc,  Tave,  overall_diff  -   100,  Tmax  =  0,  Tbv_min  =  100000,  data[N+1] [6] ; 

double  A1,  B1,  C1,  D1,  Ttip,  temp,  qa,  segment_di f f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 
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length  =  2  *  Ap  /  (  base_width  +  edge_width  ); 
b1  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2*b1  ); 
db  -  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overal l_diff )  >  (  Tbase  *  error)   | j  i  <=  N  K 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  edge_width  *  epsilon  *  pow(Ta,  4)  *  depth; 

a  =  b1  -  length; 

while  (i  <=  N  U   overflow  ==  No){ 

segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  b1  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No){ 

ua  =  u(a)  ;  ub  -   u(a  +  db); 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 

temp  =  A  *  Ta  +  B  *  qa; 

qb   =  C  *  Ta  +  D  *  qa; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  j  j  Tb  <  0){ 

overflow  -   Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

i  f    (overf low  ==   No){ 

data  [i  ][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  A; 
data[i][3]  =  B;  data[i][4]  =  C;  data[i][5]  =  D; 
Ta   =   Tb;    qa   =   qb;    a   =   a   +   db;    ++i;}} 

if    (number_of_tries   >   29){ 

clrscr();    cprintf ("Fatal    Error:      Unable   to  Compute"); 
getch();    exit(0);} 

if    (overflow  ==   NoK 

if    (Tmax   <   Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf  low  ==  YesX 

inc  =  inc  /  2; 

if  (overal l_dif f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal l_di ff; 
number_of _t  r i  es++ ; 
gotoxy(27,24); 

cprintf ("N  -   %d  Ta  -   %4 . 1 f  Tb  -   %4 . 1 f  " , number_of_tries, w_T(T , Tt ip) , w_T(T , Tb) ); 
if  (Ttip  <  Tave)  Ttip  -   Tave;}} 
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for  (i  =  N  -  1;   i  >=  1;   i--)< 

A1  =  A  *  data[i][2]  +  B  *  data[i][4];  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
C1  =  C  *  data[i][2]  +  D  *  data [i]  [4];  D1  =  C  *  data [i]  [3]  +  D  *  data[i][5]; 
A   =   A1;    B   =   B1;    C=C1;    D=D1;> 

if  (S  ==  1){ 

out  =  fopen(f i le_name,  "w"); 
for  (i=1;  i  <=  N;  i++) 

fprintf(out,  "%e,  %e, 
datati]  [0])  ); 
fclose(out);) 


%e\n",  w_L(d,  i  *  length  /  N),  w_T(T,data[i] [1] ),  w_H(H, 


Ta  =  Ttip; 

Ap  =  .5  *  length  *  (  edge_width  +  base_width  ); 

mass  =  Ap  *  depth  *  density  *  1000; 

PN  =  2  *  sigma  *  epsilon  *  pow( length,  2)  *  pow(Tbase,3)  /  (k*  base_width); 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

return;} 


/ 


*   ********************************************************** 


OUTPUT  PROCEDURE 


************************************************************************** 


*/ 


void  resul  ts(  voidH 


clrscr 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
getch( 
return 


\a\r\n 
\r\n 

\r\n 


OPTIMUM  PARAMETERS"); 
THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE"); 


Density  =  %4 

Fin  Depth  =  ") 

Total  Mass  =  ") 

Fin  Length  =  ") 

Emissivity  =  %4 

Prof i le  Area  =  ") 

Prof  i le  Number  =  %4 

Base  Heat  Flow  =  ") 

Fin  Efficiency  =  %4 

Fin  Edge  Width  =  ") 

Fin  Base  Width  =  ") 

Base  Temperature  =  ") 

Edge  Temperature  =  ") 

Thermal  Conductivity  =  ") 

Environmental  Parameter  =  %4 


,2f  kg  /  drrT3\r\n", density); 
1  p_L(L,  depth); 

p_M(M,  mass); 
•  p_L(L,  length); 
,2f \r\n",  epsi Ion); 
:  p_A(X,Ap); 
.3f\r\n",PN); 
:  p_H(H,  qb); 
,3f\r\n",eta); 

■  p_L(d,  edge_width); 

■  p_L(d,  base_width); 

■  P_T(T,Tb);  " 
'  p_T(T,Ta); 

P_K(K,  k); 
,2f\r\n",  V); 


=\r\n">; 


Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %4.4f  B  =  ",  A);  p_Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  %4.4f",  D); 


=\r\n"); 


/* 


*************  +  ***-k  +  ****±*****i 


.*************•******************** 


INPUT  PROCEDURE 
**************************************************************************  *. 


void  input(void){ 

clrscr( );  _setcursortype(_NORMALCURSOR ) ; 


cprintf ( 
cprintf ( 
cprintf ( 
cprintf ( 
cprintf ( 
cprintf ( 


\r\n 
\r\n 


THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROF  I LE\r\n"); 

Fin  Depth  ");  depth  =  i_L(L); 
Heat  Dissipation  ");  heat  =  i_H(H); 

Taper  Ratio  =  ");  scanf ("%f ",&taper_rat io); 
Emissivity  =  ");  scanf ("%f",  Sepsilon); 


=\r\n»); 


107 


cprintf ( 
cprintf ( 
cprintf ( 
cprintf  ( 
cprintf ( 


1  Density  (  kg  /  dm~3  )  =  ");  scanf("%f",  Sdensity); 

1  Base  Temperature  ");  Tbase  =  i_T(T); 

'   Thermal  Conductivity  ");  k  =  i_K(K); 

1       Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 

•\n===================================================================\r\n"); 


if  (  epsilon  >  1.0  \\    epsilon  <=  0){ 

cprintf ("\r\n\r\nFATAL  ERROR:   Emissivity  is  not  between  0  and  1"); 
getch();  exit(0);> 

if  (  V  <  0  ] |  V  >  0.8K 

cprintf ("\r\n\r\nFATAL  ERROR:   Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(0);} 

if  (  taper_ratio  >  .99  ]|  taper_ratio  <=  OX 

cprintf ("\r\n\r\nFATAL  ERROR:  Taper  Ratio  is  not  between  0  and  1.0"); 
getchO;  exit(0);} 

_set cursor type(_NOCURSOR); 
return;) 

/*  *********************************************************************** 
GLOBAL  FUNCTIONS 


******************* 


*******************************************  * 


double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));} 

double  m(void)( 

return(  hr()  /  (k  *  sin( f i n_taper) ) ) ; } 

double  u(double  x){ 

return(  2  *  sqrt(x  *  m( )  )  );) 

double  A_(void){ 

return(  ua  *  ( I 1(ua)*K0(ub)  +  I0(ub)*Kl (ua) )   );} 

double  B_(void){ 

return  (  ub*ub*( I0(ub)*K0(ua)- l0(ua)*K0(ub) )  /  (  2*k*b_width*depth*m( ) ));) 

double  C_(void){ 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*K1 (ua)  -  1 1 (ua)*K1 (ub) )/(ua)  );} 

double  D_(void)C 

return(  ub*(  I0(ua)*Kl (ub)+I 1 (ub)*K0(ua) )  );} 
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/*  *************************************************************************** 

OPTIMIZE  Analysis  Program 
Longitudinal  Fin  of  Triangular  Profile 
****************************************************************************  */ 

#include  <stdio.h> 
#include  <stdl ib.h> 
#include  <math.h> 
^include  <conio.h> 
#include  "units. h" 
#include  "bessel .h" 

#define  Yes  1 

#define  No  0 

#define  N  99 

#define  sigma  5. 66961 e- 8 

char  f i le_name [20] ; 

double  length,  edge_width,  base_width,  a_width,  b_width,  depth; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta,  PN,  heat; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  results(),  computet); 

FILE  *data_in,  *out; 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(K 

data_in  =  f open("uni ts.ext", "r"); 

fscanf (data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f ",&T,&L ,&M,SS,&H,&d,&K,SX,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

input ( ); 

heat  =  heat  /  (1  -  V); 

Ap      =  (1.1234  *  pow(heat/depth  ,  3))  /  (  k  *  pow(sigma*epsi Ion, 2)  *  pow(Tbase,9) ); 

base_width  -   (2.4852  *  pow(heat/depth  ,  2))  /  (  k  *  sigma*epsi lon*pow(Tbase,5)); 

compute( ); 
result s( ); 

ex  i  t ( 0 )  ; 
return  0;} 

/*  *********************************************************************** 

COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(){ 

double  A_(),  B_(),  C_(),  D_(),  m( ) ,  hr(),  u(double  x); 

double  length_inc,  bT,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax=0,  Tbv_min  =  100000,  data [N+1] [6] ; 

double  A1,  B1,  CI,  D1,  Ttip,  temp,  qa,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 

length  =  2  *  Ap  /  base_width; 
edge_width  -   base_width  /  100; 
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db  =  length  /  (  N  +  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*b1  ); 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001 ; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overall_diff )  >  (  Tbase  *  error)   j j  i  <=  N  K 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
segment_di f f  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

ub  =  u(db); 

temp  =  Ta  /  I0(ub); 

C    -   2  *  k  *  edge_width  *  depth  *  pow(m(),2)  *  I1(ub)  /  (ub  *  I0(ub)); 

qb   -   C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  [ j  Tb  <  OK 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  -    (Tbv_min  -  Tmax)  /  2; 

Ttip  =    fabs(Tbv_min   -    inc); 

if    (Ttip   <   Tave)    Ttip  =   Tave;}) 

i  f    (overflow  ==  NoK 

data[i][0]  =  qb;  data[i][1]  =  Tb;  data[i][2]  =  1; 
data[i][3]  =  0;  data [i]  [4]  =  C;  data[i][5]  =  1; 
Ta   =   Tb;    qa   =   qb;    ++i;) 

a   =   b1    -    length; 

while  (i  <=  N  &&  overflow  --   NoK 

segment_diff  =  10; 

a_width  -   edge_width  *  a  /  (  b1  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  bl  -  length); 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  -   A_();  B  =  B_();  C  =  C_();  D  =  D_(); 

temp  =A*Ta+B*qa; 

qb   =  C  *  Ta  +  D  *  qa; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +  50  j  J  Tb  <  OK 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbvjnin  )  Tbvjnin  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =    fabs(Tbv_min   -    inc); 

if    (Ttip  <   Tave)    Ttip  =   Tave;)) 

if    (overflow  ==   NoK 

data[i][0]  =  qb;  data[i]  [1]  =  Tb;  data[i]  [2]  =  A; 
data[i][3]  =  B;  data[i][4]  =  C;  data[i][5]  =  D; 
Ta  =  Tb;   qa  =  qb;   a  =  a  +  db;   ++i;» 
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if  (number_of_tries  >  29 ){ 

clrscrO;  cprintf ("Fatal  Error:  Unable  to  Compute"); 
getch();  exit(O);} 

if  (overflow  ==  NoH 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  YesX 

inc  =  inc  /  2; 

if  (overal l_dif f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal l_diff; 
number_of _t r i  es++ ; 
gotoxy(27,24); 

cprintf("N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tries,w_T(T, Tt ip) ,w_T(T,Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

for    (i    =   N    -    1;    i    >=    1;    )--){ 

A1  =  A  *  data[i][2]  +  B  *  data[i][A];  B1  =  A  *  data[i][3]  +  B  *  data[i][5]; 
C1  =  C  *  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i][5]; 
A   =   A1;    B   =   B1;    C  =  C 1 ;    D=D1;) 

length   =    length   +   db; 

if    (S   ==    1){ 

out    =    f open(f i le_name,    "w"); 
for    ( i    =    1;    i    <=   N;    i++) 

fprintf(out,    "%e,      %e,      %e\n",    w_L(d,    i    *    length   /   N),    w_T(T,data [i] [1] ) ,    w_H(H, 
datati]  [0]  )); 
f c  lose(out ); } 

Ta  =  Ttip; 

mass  =  Ap  *  depth  *  density  *  1000; 

PN  =  2  *  sigma  *  epsilon  *  pow( length, 2)  *  pow(Tbase,3)  /  (k*base_width); 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4) ); 

return;} 


/*  * 


********************************************************************** 


OUTPUT  PROCEDURE 


**************************************************************************  *, 


void  resul ts( voidK 


clrscr 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 
cprint 


\a 

\r\n 

\r\n 


OPTIMUM  PARAMETERS"); 
THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE"); 


:\r\n"); 


Densi  ty 

Fin  Depth 

Total  Mass 

Fin  Length 

Emi  ssi vi ty 

Prof  i le  Area 

Prof i le  Number 

Base  Heat  Flow 

Fin  Efficiency 

Fin  Base  Width 

Base  Temperature 

Edge  Temperature 

Thermal  Conductivity 

Environmental  Parameter 


=  %4.2f  kg  /  dm"3\r\n",density); 
=  ");  p_L(L,  depth); 

-  ");  p_M(M,  mass); 

-  ");  p~L(L(  length); 

=  %4.2f\r\n",  epsilon); 

=  ");  p_A(X,Ap); 

=  %4.3f\r\n",PN); 

=  ");  p_H(H,  qb); 

=  %4.3f\r\n",eta); 

=  ");  p_L(d,  base_width); 

=  ");  p_T(T,Tb); 

=  ">;  p_T(T,Ta); 

=  ");  p_K(K,  k); 

=  %4.2f\r\n",  V); 


Ar\n"); 
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cprintf (" 
cprintf (" 
cprintf(" 
cprintf (" 
getch(); 
return;} 


Transmission  Parameter  Matrix  \r\n\r\n"); 
A  =  %4.4f  B  =  ",  A);  p_Z(Z,B); 

C  =  ");  p_Y(Y,C); 
D  =  %A  . 4 f  "  ,  D); 


/*  ************************************* ********************************** 


INPUT  PROCEDURE 


**************************************************************************  *  , 


void  input (void){ 


clrscr( 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 
cprintf 


;_setcursortype(_NORMALCURSOR); 

•\r\n         THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE"); 

■\r\n    =====z============================================================\r\n"); 

1  Fin  Depth  " ) ;  depth  =  i_L(L); 

1  Heat  Dissipation  ");  heat  =  i_H(H); 

1  Emissivity  =  ");  scanf("%f",  Sepsilon); 

1  Density  (  kg  /  dnT3  )  =  ");  scanf("%f",  fcdensity); 

1  Base  Temperature  ");  Tbase  =  i_T(T); 

1    Thermal  Conductivity  ");  k  =  i_K(K); 

1        Environmental  Parameter  (0  to  .8)  =  ");  scanf("%f",  &V); 

.\n=============^^=======^===========^========^======^===r===========\r\n"); 


if  (  epsilon  >  1.0  |j  epsilon  <=  OX 

cprintf ("\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getch();  exit(O);) 

if  (  V  <  0  | |  V  >  0.8){ 

cprintf ("\r\n\r\nFATAL  ERROR:   Environmental  Parameter  is  not  between  0  and  0.8"); 
getch();  exit(0);} 

_setcursortype(_NOCURSOR ) ; 
return; ) 


****************************^ 


t************* 


GLOBAL  FUNCTIONS 
**************************************************************************  *. 

double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

return(sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));) 

double  m(void){ 

return(  hr()  /  (k  *  sin( f in_taper) ) ); ) 

double  u(double  x){ 

return(  2  *  sqrt(x  *  m( )  )  );} 

double  A_(void)( 

return(  ua  *  ( I 1(ua)*K0(ub)  +  I0(ub)*K1(ua) )   );) 

double  B_(void)<: 

return  (  ub*ub*( I0(ub)*K0(ua)- I0(ua)*K0(ub) )  /  (  2*k*b_width*depth*m( ) ) );) 

double  C_(void){ 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*Kl (ua)  -  11 (ua)*Kl (ub) )/(ua)  );) 

double  D_(void){ 

return(  ub*(  I0(ua)*Kl (ub)+I 1 (ub)*K0(ua) )  );) 
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/*  A************************************************************************** 

UNIT  CONVERSION  Program  Header  File 
A***************************************************************************  * 

#include  <stdio.h> 
#include  <conio.h> 
#inctude  <string.h> 
#include  <stdl i b. h> 
#include  "units. h" 


float  i_T(int  T){ 
float  V; 
switch  (T){ 

case  0 

case  1 

case  2 

case  3 
return( V); } 


cprintf("(deg  K)  =  ");  scanf("%f",  &V);  break; 

cprintf("(deg  C)  =  ");  scanf("%f",  &V);  V  =  V  +  273.15;  break; 

cprintf("(deg  F)  =  ");     scanf("%f",  &V);  V  =  5  *  (  V  +  459.4  )  /  9;  break; 

cprintf("(deg  R)  =  ");  scanf("%f",  &V);  V  =  5  *  V  /  9;) 


void  p_T(int  T,  float  V)( 


switch  (T)( 
case  0 
case  1 
case  2 
case  3 

return;} 


cprintf ("%4.2f  deg  K\r\n",  V);  break; 

V  =  V  -  273.15;  cpr intf ("%4 .2f  deg  C\r\n",  V);  break; 

V  =  9  *  (  V  -  255.222222)  /  5;  cprintf ("%4.2f  deg  F\r\n",  V);  break; 

V  =  9  *  V  /  5;  cprintf ("7o4.2f  deg  R\r\n",  V);> 


float  w_T(int 
switch  (TK 
case  0 
case  1 
case  2 
case  3 
return(V); ) 


float  V){ 

break; 

V  =  V  -  273.15;  break; 

V  =  9  *  (  V  -  255.222222)  /  5;  break; 
V=9*V/5;} 


void  p_Z(int  Z,  float  V)( 
switch  (Z){ 

case  0:      cprintf ("%4.4f  deg  C   /  watt\r\n",  V);  break; 

case  1:      V  =  0.1760831  *  V  ;  cpr intf ("%4.34  deg  F  hr  /  BTU\r\n",  V);} 
return;) 

float  w_Z(int  Z,  float  V){ 
switch  (Z)( 

case  0:      break; 

case  1:      V  =  0.1760831  *  V;) 
return( V); ) 

void  p_Y(int  Y,  float  V)( 
switch  (Y){ 

case  0:      cprintf ("%/.  .4f  watt   /  deg  C",  V);  break; 

case  1:      V  =  V  /  0.1760831  ;  cprintf ("%A.4f  BTU  /  hr  deg  F",  V);} 
return;} 

float  w_Y(int  Y,  float  V){ 
switch  (Y){ 

case  0: 

case  1 : 
return(V); } 

float  i_L(int  L){ 
float  V; 
switch  (L)( 

case  0 

case  1 

case  2 

case  3 

case  A 


break; 

V  =  V  /  0.1760831;} 


cprintfC"  (m)  =  ");  scanf("7of 

cprintf("(cm)  =  ");  scanf("7.f 

cprintf("(ft)  -   ");  scanf("%f 

cprintf  ("(in)  -   ");  scanf("7.f 

cprintf ("(mm)  -   ");  scanf("%f 


".    &V); 

",  &v); 

",    &V); 

",    &V); 

",    &V); 

break; 

V  =  V  /  100;  break; 

V  =  0.3048  *  V;  break; 

V  =  0.0254  *  V;  break; 

V  =  V  /  1000;} 
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return(V);} 


void  p_L( int  L 
switch  (L){ 
case  0 
case  1 
case  2 
case  3 
case  4 
return;} 


float  V)( 


cprrntf("%4.2f  m\r\n",  V);  break; 

V  =  100  *  V;  cprintf("%4.2f  cm\r\n",  V);  break; 

V  =  V  /  0.3048;  cprintf ("%4.2f  ft\r\n",  V);  break; 

V  =  V  /  0.0254;  cprintf ("%4.2f  in\r\n",  V);  break; 

V  =  1000  *  V;  cprintf("%4.2f  rmt\r\n",  V);) 


float  w_L(int  L,  float  V)C 


switch  (LH 
case  0 
case  1 
case  2 
case  3 
case  4 

return(V);} 


break; 

V  =  100  *  V;  break; 

V  =  V  /  0.3048;  break; 

V  =  V  /  0.0254;  break; 

V  =  1000  *  V;> 


float  i_A(int  A){ 
float  V; 
switch  (A){ 

case  0 

case  1 

case  2 

case  3 

case  4 
return(V);} 


cprintfC  (m"2) 
cprintf  C(cnT2) 
cprintf("(ft"2) 
cprintf ("( in"2) 
cprintf ("(mm"2) 


=  ");  scanf("%f", 

=  ");  scanf("%f", 

=  ");  scanf("%f", 

=  ");  scanf("%f", 

=  ");  scanf("7.f", 


&V); 
&V); 


break; 

V  =  V  / 
I 


10000;  break; 
10.76;  break; 

V  =  V  /  1550;  break; 

V  =  V  /  1000000;) 


&V);  V  =  V 

&V); 

&V); 


void  p_A(int  A,  float  V){ 


switch  (A){ 
case  0 
case  1 
case  2 
case  3 
case  4 

return; } 


cprintf ("%4. 2  m  "  2\r\n",  V);  break; 

V  =  10000  *  V;  cprintf("%4.2f  cm  "  2\r\n",  V);  break; 

V  -   V  *  10.76;  cprintf("%4.2f  ft  '  2\r\n",  V);  break; 

V  -   V  *  1550;  cprintf("%4.2f  in  "  2\r\n",  V);  break; 

V  =  1000000  *  V;  cprintf("7.4.2f  nvn  '  2\r\n",  V);) 


float  w_A(int  A,  float  V){ 


switch  (A){ 
case  0 
case  1 
case  2 
case  3 
case  4 

return(V);} 


break; 

V  -    10000  *  V;  break; 

V  =  V  *  10.76;  break; 

V  =  V  *  1550;  break; 

V  =  1000000  *  V;} 


float  i_M(int  M){ 
float  V; 
switch  (M){ 

case  0 

case  1 

case  2 

case  3 
return(V);} 


cprintfC  (kg)  =  ");  scanf("%f",  &V);  break; 

cprintfC  (gm)  =  ");  scanf("%f",  &V);  V  =  V  /  1000;  break; 

cprintf C(lbm)  =  ");  scanf("%f",  &V);  V  =  V  /  2.204622;  break; 

cprintfC  (oz)  =  ");  scanf ("%f ",  &V);  V  =  V  /  35.27396;  break;} 


void  p_M( int  M,  float  V){ 
switch  (M){ 

case  0 

case  1 

case  2 

case  3 
return;) 


cprintf("%4.2f  kg\r\n",  V);  break; 

V  =  1000  *  V;  cprintf("%4.2f  grams\r\n",  V);  break; 

V  -   2.204622  *  V;  cprintf ("%4.2f  lbm\r\n",  V);  break; 

V  =  35.27396  *  V;  cprintf ("%4.2f  oz\r\n",  V);} 


float  w_M(int  M,  float  V)( 
switch  (M){ 

case  0:      break; 

case  1:      V  =  1000  *  V;  break; 
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case  2: 
case  3: 
return(V); } 

float  i_H(int  H)( 
float  V; 
switch  (H){ 

case  0 

case  1 

case  2 

case  3 
return(V);) 


V  =  2.204622  *  V;  break; 

V  =  35.27396  *  V;) 


cprintf("(watts)  =  ");  scanf("%f",  &V);  break; 
cprintf("     (kw)  =  ");  scanf("%f",  &V);  V  =  1000  *  V;  break; 
cprintf("(BTU/sec)  =  ");  scanf("%f",  &V);  V  =  V  *  1054.3502;  break; 
cprintf("  (BTU/hr)  =  ");  scanf("%f",  &V);  V  =  V  /  3.414425;) 


void  p_H(int  H,  float  V){ 
switch  (H){ 

case  0 

case  1 

case  2 

case  3 
return;} 


cprintf("%4.2f  watts\r\n",  V);  break; 

V  =  V  /  1000;  cprintf("%4.2f  kw\r\n",  V);  break; 

V  =  V  /  1054.3502;  cprintf ("%4.2f  BTU  /  sec\r\n",  V);break; 

V  =  3.414425  *  V;  cprintf ("%4.2f  BTU  /  hr\r\n",  V);) 


float  w_H(int  H,  float  VK 


switch  (H)C 
case  0 
case  1 
case  2 
case  3 

return(V);} 


break; 

V  =  V  /  1000;  break; 

V  =  V  /  1054.3502;  break; 

V  =  3.414425  *  V;} 


float  i_K< int  K)( 

float  V; 

switch  (K){ 

case  0: 
case  1 : 

return(V);} 


cprintf ("(watts  /  m  deg  K)  =  ");  scanf("%f",  &V);  break; 
cprintf("(BTU  /  ft  hr  deg  F)  =  ");  scanf("%f",  &V);  V  =  V  /  1.731;} 


void  p_K(int  K,  float  V){ 
switch~(K)( 

case  0:      cprintf ("%4.2f  watts  /  m  deg  K\r\n",  V);  break; 

case  1:      V  =  1.731  *  V;  cprintf ("%4.2f  BTU  /  ft  hr  deg  F\r\n",  V);} 
return; } 
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/*  ************************************************************************** 

BESSEL  FUNCTION  Program  Header  File"' 
/*  *************************************************************************** 

#include  <stdio.h> 
#include  <math.h> 
#inc lude  "bessel .h" 

double  K0(float  x){ 

double       y,  z; 

double       Pl=-0. 57721566,  P2=0. 4227842,  P3=0. 23069756,  P4=0. 0348859,  P5=0. 00262698, 

P6=0. 0001075,  P7=0. 0000074; 
double       01=1.25331414,  Q2=-0. 07832358,  03=0.02189568,  Q4=-0. 01062446,  05=0.00587872, 

Q6=-0. 00251546; 
double       07=0.00053208; 

if  (  x  <  2  ){ 

y  =  (  x  *  x  )  /  4; 

z  =  (-log(x/2)*I0(x))+(P1+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7) )))));) 
else{ 

y  =  2  /  x; 

z  =  (exp(-x)/sqrt(x))*(Q1+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*Q7))))));> 
return( z) ; } 

double  KKfloat  x){ 

double       y,  z; 

double       P1=1.0,  P2=0. 15443144,  P3=- .67278579,  P4=-0. 18156897,  P5=-0. 01919402,  P6=-0. 001 10404, 

P7=-0. 00004686; 
double       01=1.25331414,  02=0.23498619,  Q3=-0. 03655620,  04=0.01504268,  Q5=-0. 00780353, 

06=0.00325614; 
double       Q7=-0. 00068245; 

if  (  x  <  2  X 

y  =  (  x  *  x  )  /  4; 

z  =  (log(x/2)*H(x))+(1/x)*(P1+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7))))));> 
else{ 

y  =  2  /  x; 

z  =  (exp(-x)/sqrt(x))*(Ql+y*(Q2+y*(Q3+y*(04+y*(Q5+y*(06+y*Q7) )))));) 
return(z);) 

double  I0(float  x){ 

double       y,  z,  ax; 

double       P1=1,  P2=3. 5156229,  P3=3. 0899424,  P4=1 .2067492,  P5=0. 2659732,  P6=0. 0360768, 

P7=0. 0045813; 
double       01=0.39894228,  02=0.01328592,  03=0.00225319,  Q4=-0. 00157565,  05=0.00916281, 

06=- 0. 02057706; 
double       07=0.02635537,  08= -0. 01647633,  09=0.00392377; 

if  (fabs(x)  <  3.75  )( 

y  =  (  x  /  3.75  )  *  (  x  /  3.75  ); 

z  =  P1+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7)))));) 
else{ 

ax  =  fabs(x); 

y  =  3.75  /  ax; 

z  =  (exp( ax )/sqrt(ax))*(Q1+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*(Q7+y*(Q8+y*Q9) )))))));> 
return(z);} 

double  IKfloat  x){ 

double       y,  z,  ax; 

double       P1=0.5,  P2=0. 87890594,  P3=0. 51498869,  P4=0. 15084934,  P5=0. 02658733,  P6=0. 00301532, 

P7=0. 000324 11; 
double       01=0.39894228,  Q2=-0. 03988024,  Q3=-0. 00362018,  04=0.00163801,  Q5=-0. 01031555, 
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06=0.02282967; 
double       Q7=-0. 02895312,  08=0.01787654,  Q9=-0. 00420059; 

if  (fabs(x)  <  3.75  )( 

y  =  (  x  /  3.75  )  *  (  x  /  3.75  ); 

z  =  x*(PHy*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7))))));) 
else{ 

ax  =  fabs(x); 

y  =  3.75  /  ax; 

Z  =  (exp(ax)/sqrt (ax) )*(Ql+y*(02+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*(Q7+y*(Q8+y*Q9) )))))));> 
return(z);} 
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APPENDIX  B  COMPUTER  CODE  TO  DERIVE  OPTIMUM  GEOMETRY 

/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Rectangular  Profile 
****************************************************************************  */ 

#include  <time.h> 
# include  <stdio.h> 
tfinclude  <stdl ib.h> 
#include  <math.h> 
#include  <conio.h> 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

double  length,  width_base,  max_length; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  Heat,  qb_max,  db,  Continue; 

float  V,  K2,  epsilon,  k,  Tbase,  Ap,  taper_ratio,  depth,  error=0.01; 

int  i,  j,  number_of_passes,  failures; 

void  computet ); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(){ 

float  dq; 

FILE  *out; 

randomi  ze( ); 

error  =  error  /  100; 

clrscr( ); 

for  (j=1;  j<102;  j++){ 

number_of_passes  =  0; 

epsilon  =  (float)  (  random(98)  ) ; 

epsilon  =  .01  +  epsilon  /  100; 

depth  =  (float)  (20  +  random(200)); 

depth  =  depth  /  100; 

k  =  10  +  random(200); 

Tbase  =  100  +  random(500); 

Heat  =  (float)  (25  +  random(IOOO)); 

Heat  =  Heat  /  depth; 

V  =   (float)   (    random(80)); 

V  =    .01    +   V   /    100; 

V  =   0; 

taper_ratio  =    (float)    (random(98) ); 
taper_ratio  =   0.01    +   taper_ratio  /    100; 
taper_ratio   =    1; 

Continue  =   Yes; 
qb_max   =    -1; 
qb  =   0; 
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dq  =  0; 

Ap  =  1.635  *  pow(Heat,3)  /  (  k  *  pow(sigma,2)  *  pow(epsi Ion, 2)  *  pow(Tbase,9) ); 

length  =  .87  *  Heat  /  (  sigma  *  epsilon  *  pow(Tbase,4)); 

db  =  length  /  20016; 

length  =  length  -  db; 

width_base  =  Ap  /  length; 

if  (depth  <  10*width_base)  Continue  =  No; 
if  (length  <  10*width_base)  Continue  =  No; 

while  (Continue  ==  Yes)( 

length  =  length  +  db; 
width_base  =  Ap  /  length; 
number_of_passes++; 
fai lures  =  0; 

compute( ); 

while  (qb  ==  0  &&  Continue  ==  Yes){ 
fai lures++; 

length  =  length  +  db  /  2; 
width_base  =  Ap  /  length; 
c  ompu  t  e  ( ) ; 
if  (fai lures  ==  3)( 

Continue  =  No; 

number_of_passes  =  1;}> 

dq  =  qb  -  dq; 

cprintf("\r\nwidth_base  =  %A.3f  length  =  %A.6f  heat  =  %A.6f  dq  =  %4.6f", 

width_base*1000,  length,  qb,  dq); 
dq  =  qb; 

if  (qb_max  <  qb  ){ 

max_length  =  length; 

qb_max  =  qb;) 
else 

Continue  =  No;) 

if  (number_of_passes  >  2){ 

length  =  max_length; 

width_base  =  Ap  /  length; 

compute( ); 

cprintf("\r\nWriting  Rec  #  %d  d  =  %A.3f  b  =  %4.6f  qb  =  %4.3f  Tb  =  %A.1f  Ap  =  %4.1f", 
j,  width_base,  length,  qb,  Tbase,  Ap); 

out  -   fopen("OPT_RECT.MAT","a"); 

fprintf(out,  "%d,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e\r\n",  j, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_ratio,  width_base,  width_base  ), 

fclose(out );) 
else{ 

cprintf ("\r\nIgnoring  This  Set  Of  Parameters"); 

J=j-1;)> 

return  0;) 

/*  a********************************************************************** 

COMPUTE  ROUTINE 
**************************************************************************  */ 

void  compute(void){ 

double       Tav,  inc,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000; 
double      A,  B,  C,  D,  m,  hr,  Yo,  Ttip,  temp,  qa,  Tave,  segment_dif f ; 
int         overflow,  possible_overf low  =  No,  number_of_tries=0; 

width_base  =  Ap  /  length; 
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K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4) ; 
Tave  -  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2), 
else 

Ttip  =  .0001  +  Tave; 

inc  =  .5  *  Ttip;  i  =  1; 

while  (fabs(overal l_dif f )  >  (  Tbase  *  error)   |]  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  -   sigma  *  width_base  *  epsilon  *  pow(Ta,  4)  *  depth; 

while  (i  <=  N  &&  overflow  ==  NoK 

segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoK 

Tav  =  (Tb  +  Ta)  /  2; 

hr  =  sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)); 

m  =  sqrt(  (  2.0  *  hr)  /  (  k  *  width_base)  ); 

Yo  =  sqrt(  2.0  *  hr  *  k  *  width_base)  *  depth; 

A  =  cosh(  m  *  length  /  N  ); 

B  =  sinh(  m  *  length  /  N  )  /  Yo; 

C  =  sinh(  m  *  length  /  N  )  *  Yo; 

D  =  A; 

temp  =A*Ta+B*qa; 
qb   =  C  *  Ta  +  D  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  J  j  Tb  <  0){ 

overflow  -   Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;)} 

if  (overflow  ==  NoK 

Ta  =  Tb;  qa  =  qb;  ++i;)} 

if  (number_of_tries  >  29){ 

Ta=0;Tb=0;qb=0; return;} 

if  (overflow  ==  No)( 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase  )  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  --   Yes){ 
inc  =  inc  /  2 ; 

if  (overal L_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  +  inc; } 
else 

Ttip  =  Ttip  -  overal l_di ff; 
number_of _t r i  es++  ; 
if  (Ttip  <  Tave)  Ttip  =  Tave;)} 
return;) 
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/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program  To  Determine  The  Optimum  Geometry 

Longitudinal  Fin  of  Trapezoidal  Profile 
****************************************************************************  */ 

#include  <time.h> 
#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
^include  <conio.h> 
#include  "bessel.h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

double  max_length,  length,  depth,  edge_width,  base_width,  Ta,  Tb,  qb=0,  inc,  Tave; 

float  dq,  V,  qbjnax,  db,  epsilon,  taper_ratio,  k,  error=0.01,  Tbase,  K2,  Heat,  Ap; 

int  i,  Continue,  number_of_passes,  j,  Exit,  failures; 

void  computet); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

i  nt  mai  n( )( 

FILE  *out; 
randomi  ze(  ); 
error  =  error  /  100; 
clrscr( ); 

for  (j=1;  j  <  102;  j++H 

number_of_passes  =  0; 

epsilon  =  (float)  (  random(98)  ); 

epsilon  =  .01  +  epsilon  /  100; 

depth  =  (float)  (20  +  random(200) ); 

depth  =  depth  /  100; 

k  =  10  +  random(200); 

Tbase  =  100  +  random(500); 

Heat  =  (float)  (25  +  random(IOOO)); 

Heat  =  Heat  /  depth; 

V  =  (float)  (  random(80)); 

V  -  .01  +  V  /  100; 

V  =  0; 

taper_ratio  =  (float)  (random(98) ); 
taper_ratio  =  0.01  +  taper_ratio  /  100; 

Continue  =  Yes; 

Exit  =  No; 

qb_max  =  - 1 ; 

dq  =  0; 

qb  =  0; 

Ap  =  (.63077  *  taper_ratio  +  1.31723)  *  pow(Heat,3)  /  (  k  *  pow(sigma,2)  *  pow(epsi Ion, 2) 

pow(Tbase,9)); 
length  =  .82  *  (-.0315113  *  taper_ratio  +  1.0425113)  *  Heat  /  (  sigma  *  epsilon  * 

pow(Tbase,4) ); 
db  =  length  /  20016; 

length  =  length  -  db; 

base_width  =  2  *  Ap  /  (  length  *  (  taper_ratio  +  1)); 

edge_width  =  taper_ratio  *  base_width; 
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if  (depth  <  10*base_width)  Continue  =  No; 
if  (length  <  10*base_width)  Continue  =  No; 

while  (Continue  ==  Yes)( 

length  =  length  +  db; 

base_width  -   2   *  Ap  /  (  length  *  (  taper_ratio  +  1)); 

edge_width  =  taper_ratio  *  base_width; 

fai lures  =  0; 
number_of_passes++; 

compute(  ); 

while  (qb  ==  0  &&  Continue  ==  YesH 
fai lures++; 

length  =  length  +  db  /  2; 

base_width  =  2  *  Ap  /  (  length  *  (  taper_ratio  +  1)); 
edge_width  =  taper_ratio  *  base_width; 
compute ( ); 
if  (failures  ==  3){ 

Continue  -   No; 

number_of_passes  =  1;)> 

dq  =  qb  -  dq; 

cprintf("\r\nwidth  =  %4.3f  length  =  %A.6f  heat  =  %A.6f  dq  =  %A.6f",  base_width*1000, 

length,  qb,  dq); 
dq  =  qb; 

if  (qb_max  <  qb  ){ 

max_length  =  length; 

qb_max  =  qb;} 
else 

Continue  =  No;  > 

if  (number_of_passes  >  2)( 

length  =  max_length; 

base_width  =  2  *  Ap  /  (  length  *  (  taper_ratio  +  1)); 

edge_width  =  taper_ratio  *  base_width; 

compute ( ) ; 

cprintf ("\r\nWriting  Rec  #  %d  d  =  %4.3f   b  =  %4.6f  qb  =  %A.3f  Tb  -   %4.1f  Ap  =  %4.1f", 
j,  base_width,  length,  qb,  Tbase,  Ap); 

out  =  fopen("OPTTRAP.MAT","a"); 

fprintf(out,  "%d,   %e,  %e,  %e,  %e,  %e,  %e,  %e,   %e,   %e,  %e,   %e\r\n",  j, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_ratio,  base_width,  edge_width  ); 

f close(out ); } 
else{ 

cprintf ("\r\nlgnoring  This  Set  Of  Parameters"); 

J=j-1;)> 

return  0;> 

void  compute(void){ 

double  m,  hr,  n,  ua,  ub,  A,  B,  C,  D,  fin_taper,  I1UA,  KOUB,  I0UB,  K1UA,  KOUA,  IOUA,  K1UB; 

double  qa,  b1 ,  dH,  a,  a_width,  b_width,  Ttip,  temp,  segment_di f f ,  HUB; 

double  b,  inc,  Tav,  Tave,  overall_diff  -    100,  Tmax  =  0,  Tbvjnin  =  100000; 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 

b1  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2  *  b1  ); 
dH  =  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4) ; 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
else 
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Ttip  =  Tave  +  0.0001; 
inc  =  .5  *  Ttip;  i=1; 
while  (f abs(overal I  diff)  >  (  Tbase  *  error) 


i  <=  N  ){ 


overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  edge_width  *  epsilon  *  pow(Ta,  4)  *  depth; 

a  =  b1  -  length; 

while  (i  <=  N  &&  overflow  —  NoX 

segment_dif f  =  10; 

a_width  a  edge_width  *  a  /  (  bl  -  length); 

b_width  -   edge_width  *  (  a  +  dH  )  /  (  bl  -  length); 

while  (fabs(segment_di f f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  NoX 

Tav  =  (Tb  +  Ta)  /  2; 

hr  =  sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)); 

n  =  hr  /  (k  *  sin(f in_taper)); 

ua  -     2  *  sqrt(a  *  n  ); 

ub  -      2  *  sqrt(  (a  +  dH)  *  n  ); 

I1UA  =  IKua); 

K0UB  -  K0(ub); 

I0UB  -    I0(ub); 

K1UA  =  Kl(ua); 

K0UA  =  K0(ua); 

I0UA  =  I0(ua); 

K1UB  -   KKub); 

HUB  =  IKub); 

A  -  ua  *  (I1UA*K0UB  +  10UB*K1UA); 

B  =  ub*ub*(10UB*K0UA  -  I0UA*K0UB)  /  (  2*k*b_width*depth*n); 

C  =  2*k*a_width*depth*n*ub*(HUB*KlUA  -  I lUA*KlUB)/(ua); 

D  =  ub*(  !0UA*K1UB  +  11UB*K0UA); 
temp  =  A  *  Ta  +  B  *  qa; 
qb   =  C  *  Ta  +  D  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  j  ]  Tb  <  OX 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

if  (overflow  ==  NoX 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  dH;  ++i;}} 

if  (number_of_t r les  >  29)( 

Ta=0;Tb=0;qb=0; return;) 

if  (overflow  ==  NoX 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf  low  ==  YesX 

inc  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;} 
else 

Ttip  =  Ttip  -  overal l_di ff; 
number_of _t  r i  es++ ; 
if  (Ttip  <  Tave)  Ttip  =  Tave;}} 


return;} 
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/*  *************************************************************************** 

PROFILE  AREA  Analysis  Program  Tb  Determine  The  Optimum  Geometry 

Longitudinal  Fin  of  Triangular  Profile 
****************************************************************************  */ 

#include  <time.h> 
# include  <stdio.h> 
#include  <stdl ib.h> 
# include  <math.h> 
#include  <conio.h> 
#include  "bessel.h" 

#define  Yes  1 

#define  No  0 

tfdefine  N  99 

#define  sigma  5. 66961 e- 8 

double  length,  depth,  edge_width,  base_width,  max_length; 

double  Ta,  Tb,  qb=0; 

float  V,  qb_max,  db,  epsilon,  k,  error=0.01,  taper_ratio,  Tbase,  K2,  Heat,  Ap; 

int  i,  Exit,  Continue,  number_of_passes,  j,  failures; 

void  computet); 

/*  *********************************************************************** 

MAIN  PROGRAM 
**************************************************************************  */ 

int  main(){ 
float  dq; 

FILE  *out; 
randomi  ze( ); 
error  =  error  /  100; 
clrscr( ); 

for  (j=1;  j  <  102;  j++X 

number_of_passes  =  0; 

epsilon  -    (float)  (  random(98)  ); 

epsilon  =  .01  +  epsilon  /  100; 

depth  =  (float)  (20  +  random(200) ); 

depth  =  depth  /  100; 

k  =  10  +  random(200); 

Tbase  =  100  +  random(500); 

Heat  =  (float)  (25  +  random(IOOO)); 

Heat  -   Heat  /  depth; 

V  =  (float)  (  random(80)); 

V  =  .01  +  V  /  100; 

V  =  0; 

taper_ratio  =  (float)  (random(98)); 
taper_ratio  -   0.01  +  taper_ratio  /  100; 
taper_ratio  =  0; 

Exit  =  No; 

Continue  =  Yes; 

qb_max  =  - 1 ; 

qb  =  0; 

dq  =  0; 

Ap  =  1.304  *  pow(Heat,3)  /  (  k  *  pow(sigma,2)  *  pow(epsi Ion, 2)  *  pow(Tbase,9) ); 

length  =  .95  *  Heat  /  (  sigma  *  epsilon  *  pow(Tbase,4) ); 

db  =  length  /  20016; 

length  =  length  -  db; 

base_width  -   2  *  Ap  /  length; 
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if  (depth  <  10*base_width)  Continue  =  No; 
if  (length  <  10*base_width)  Continue  =  No; 

while  (Continue  ==  YesK 

length  =  length  +  db; 
base_width  =  2  *  Ap  /  length; 

number_of_passes++; 
fai lures  =  0; 

compute ( ); 

while  (qb  ==  0  &&  Continue  ==  YesK 
fai lures++; 

length  =  length  +  db  /  2; 
base_width  =  2  *  Ap  /  length; 
compute( ); 
if  (fai lures  ==  3)( 
Continue  =  No; 
number_of_passes  =  1;)> 

dq  =  qb  -  dq; 

cprintf("\r\ndelta  =  %4.3f  length  =  %A.6f  heat  =  %4.6f  dq  =  %4.6f",  base_width*1000, 

length,  qb,  dq); 
dq  =  qb; 

if  (qb_max  <  qb  ){ 

max_length  -  length; 

qb_max  =  qb; > 
else 

Continue  =  No;  ) 

if  (number_of_passes  >  2){ 

length  =  max_length; 

base_width  =  2  *  Ap  /  length; 

compute( ); 

cprintf("\r\nWriting  Rec  #  %d  d  =  %A.3f  b  =  %4.6f  qb  =  %4.3f  Tb  =  %4 . 1 f  Ap  =  %4.1f", 
j,  base_width,  length,  qb,  Tbase,  Ap); 

out  =  fopen("OPTTRI .MAT","a"); 

edge_width  -   0; 

fprintf(out,  "%d,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e\r\n",  j, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_ratio,  base_width,  edge_width  ); 

fclose(out);} 
else{ 

cprintf ("\r\nlgnoring  This  Set  Of  Parameters"); 

J=j-1;» 

return  0;) 

void  compute(void){ 

double  m,  hr,  A,  B,  C,  D,  I1UA,  KOUB,  IOUB,  K1UA,  KOUA,  IOUA,  K1UB,  HUB; 

double  fin_taper,  b1 ,  db,  a,  ua,  ub,  a_width,  b_width; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbvjnin  =  100000; 

double  Tav,  Ttip,  temp,  qa,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 

edge_width  =  base_width  /  100; 

db  =  length  /  (  N  +  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*b1  ); 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 
Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 
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else 

Ttip  =  Tave  +  0.0001; 

inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overal l_dif f )  >  (  Tbase  *  error)   j |  i  <=  N  )( 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  -   Ta; 
segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No)( 

Tav  =  (Tb  +  Ta)  /  2; 

hr  =  sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)); 
m  =  hr  /  (k  *  sin(f in_taper)); 
ub  =  2  *  sqrt(  db  *  m  ); 
I0UB  =  I0(ub), 
K1UB  =  KKub), 
HUB  =  IKub), 
C  =  2*k*edge_width*  depth*pow(m,2)*I lUB/(ub  *  I0UB); 


temp  =  Ta  /  I0UB; 

qb   =  C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  +50  j  j  Tb  <  0)( 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 

inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;}} 

if  (overflow  ==  No){ 

Ta  =  Tb;  qa  =  qb;  ++i ;} 

a  =  b1  -  length; 

while  (i  <=  N  &&  overflow  ==  No){ 

segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  b1  -  length); 

b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (f abs(segment_di f f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No){ 


Tav  =  (Tb  +  Ta)  /  2; 

hr  =  sigma  *  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)); 

m  =  hr  /  (k  *  si n( f i n_taper) ); 

ua  =  2  *  sqrt(a  *  m  ); 

ub  -     2  *  sqrt(  (a  +  db)  *  m  ); 

I1UA  =  IKua); 

K0UB  =  K0(ub); 

I0UB  =  10(ub); 

K1UA  =  KKua); 

K0UA  =  K0(ua); 

I0UA  =  I0(ua); 

K1UB  =  KKub); 

HUB  =  IKub); 

A  =  ua  *  (I1UA*K0UB  +  I0UB*K1UA); 

B  =  ub*ub*(I0UB*K0UA  -  I0UA*K0UB)  /  (  2*k*b_width*depth*m); 

C  =  2*k*a_width*depth*m*ub*(HUB*K1UA  -  HUA*KlUB)/(ua); 

D  =  ub*(  I0UA*K1UB  +  I1UB*K0UA); 


temp  =  A  *  Ta  +  B  *  qa; 
qb   =  C  *  Ta  +  D  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
Tb  =  temp  +  segment_dif f ; 
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if  (Tb  >  Tbase  +50  j  j  Tb  <  OX 

overflow   =  Yes;  possible_overf low  =  Yes; 

if  (Ttip  <  Tbv_rnin  )  Tbvjnin  =  Ttip; 

inc  =  (Tbvjnin  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 

if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

i  f  (overflow  ==  No){ 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  db;  ++i;}) 

if  (number_of_tries  >  29 )( 

Ta=0;  Tb=0;  qb=0;  return;} 

if  (overflow  ==  No)< 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  YesK 

inc  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  -   Ttip  -  inc; 

else  Ttip  =  Ttip  +  inc;) 
else 

Ttip  =  Ttip  -  overal l_diff; 
number_of _t r i  es++ ; 
if  (Ttip  <  Tave)  Ttip  -   Tave;)) 

length  =  length  +  db; 
return;} 
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APPENDIX  C  MATLAB  EVALUATION  M  FILES 

/*  *************************************************************************** 

MATLAB  EVALUATION  M  FILE  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Rectangular  Profile 
****************************************************************************  *  / 

load  rect; 
D  =  rect; 

N  =  100; 

sigma  =  5.66961e-8; 

L     =  D(:,2) 
b     =  D(:,3) 
A     =  D(:,4) 
k     =  D(:,5) 
To  =  D(:,6) 
e     =  D(:,7) 
V     =  D(:,8) 
qo  =  D(:,9) 
TR   =   D(:,10) 
db  =  D(:,11) 
de  =  DC : , 12) 
I      =    [1:101] 

C3   =  A.*   k.*    (((sigma*e).A2).*   To.A9).*    (L./qo).A3; 

C1    =   b.*   sigma.*   e.*    (To. '4).*    (L./qo); 

C2   =  db.*   k.*   sigma.*   e.*    (To. '5).*    (L./qo).A2; 

dispCRECT   V   =   0') 

dispC    ') 

disp( 'Constant   One:      Polynomial    p   is:') 

p  =   polyfit(TR,    C1,    1) 

fitl    =   polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum((C1-f i t1 ) . "2)/(N- 1 )) 

pause 

dispC    ') 

di sp( 'Constant   Two:      Polynomial    p   is:') 

p  =  polyfit(TR,    C2,    1) 

fit2  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum((C2-f i t2) . A2)/(N- 1 ) ) 

pause 

dispC    ') 

disp( 'Constant   Three:      Polynomial    p   is:') 

p  =   polyfit(TR,    C3,    1) 

fit3  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err   =   sqrt(sum( (C3-f i t3) . A2)/(N-1 )) 


128 


/*  ************************ *************************************************** 

MATLAB  EVALUATION  M  FILE  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Trapezoidal  Profile 
****************************************************************************  */ 

load  trap; 
D  =  trap; 

N  =  100; 

sigma  =  5.66961e-8; 

L     =  D(:,2) 
b     =  D(:,3) 
A     =  D(:,4) 
k     =  DC, 5) 
To  =  D(:,6) 
e     =  D(:,7) 
V     =  D(:,8) 
qo  =  D(:,9) 
TR   =  D(:,10) 
db  =  D(:,11) 
de  =  DC : , 12) 
I      =    [1:101] 

C3   =  A.*   k.*    (((sigma*e)."2).*   To. "9).*    (L./qo).A3; 

C1    =  b.*   sigma.*   e.*    (To. '4).*    (L./qo); 

C2   =   db.*   k.*   sigma.*   e.*    (To. "5).*    (L./qo). *2; 

dispCTRAP   V   =   0'  ) 

dispC    ') 

disp( 'Constant   One:      Polynominal    p   is:1) 

p  =  polyfit(TR,    C1,    1) 

fitl    =  polyval(p.TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum((C1 -f i t1 ) . *2)/(N- 1 ) ) 

pause 

dispC  ') 

disp( 'Constant  Two:  Polynominal  p  is:') 

p  =  polyfitCTR,  C2,  2) 

fit2  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum((C2-f i t2) . A2)/(N-1 )) 

pause 

dispC  ') 

disp( 'Constant  Three:  Polynominal  p  is:') 

p  =  polyfitCTR,  C3,  1) 

fit3  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err   =   sqrt(sum( (C3-f i t3) . *2)/(N- 1 ) ) 

pause 
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/*  *************************************************************************** 
MATLAB  EVALUATION  M  FILE  To  Determine  The  Optimum  Geometry 

Longitudinal  Fin  of  Triangular  Profile 
****************************************************************************  * 


/ 


I oad  t  r i ; 
D  =  tri; 

N  =  100; 

sigma  =  5.66961e- 


L  =  D(:,2) 
b  =  D(:,3) 
A  =  D(:,4) 
k  =  D(:,5) 
To  =  D(:,6) 
e  =  D(:,7) 
V  =  D(:,8) 
qo  =  D(:,9) 
TR  =  D(:,10) 
db  =  D(:,11) 
de  =  D(:,12) 
I   =  [1:101] 


C3  =  A.*  k.*  (((sigma*e).~2).*  To. "9).*  (L./qo).*3; 

C1  =  b.*  sigma.*  e.*  (To. "A).*  (L./qo); 

C2  =  db.*  k.*  sigma.*  e.*  (To. "5).*  (L./qo). "2; 

disp( 'TRI  V  =  0' ) 

dispC  ') 

disp( 'Constant  One:   Polynominal  p  is:') 

p  =  polyf it(TR,  C1,  1) 

fitl  -   polyval(p.TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum( (C1 -f i t 1 ) . ~2)/(N- 1 ) ) 

pause 

dispC  ') 

disp( 'Constant  Two:  Polynominal  p  is:') 

p  =  polyfit(TR,  C2,  1) 

fit2  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err  =  sqrt(sum((C2-f i t2) . A2)/(N- 1 )) 

pause 

dispC  ') 

disp( 'Constant  Three:  Polynominal  p  is:') 

p  =  polyf it(TR,  C3,  1) 

fit3  =  polyval(p,TR); 

disp('Mean  Square  Error:') 

std_err   =   sqrt(sum((C3-f i t3)."2)/(N-1 )) 
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